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result  in  practice  as  well  as  in  theory, 
practice.  _ ^ 


This  report  shows  this  equivalence  in 


JThe  atmosphere  is  assumed  to  consist  of  two  layers.  The  lower  layer, 
which  is  immediately  above  the  earth's  surface,  is  a  non-nuctiug  one,  typically 
the  standard  atmosphere,  with  a  constant  thickness,  and  the  upper  one  is  a 
ducting  layer  of  sufficient  thickness.  On  the  condition  that  the  height  above 
the  earth's  surface  and  the  elevation  angle  of  the  ray  remain  sufficiently  small 
throughout  the  ray  trajectory,  it  can  be  shown  that  the  ray  path  in  either  layer 
may  be  approximated  as  a  circular  arc  of  appropriate  radius.  \  Using  this 
condition,  one  can  derive  a  set  of  equations  that  can  be  used  to)  determine  the 
shape  of  the  ray  paths.  The  equation  for  the  angle  of  take-off pr  arrival  that 
follows  from  the  set  and  comes  in  the  form  of  a  quartic  equajfon  can  sometimes 
have  more  than  one  acceptable  solution.  The  collectionpMhese  solutions  that 
have,  in  general,  different  angles  and  times  of  arjjxrfw ith  different  relative 
intensities  can  be  considered  to  represent .thglffultipath  for  the  link  under 
consideration.  -  — -  — ~ 

^  A  comparison  of  our  quartic  equation  from  the  Earth  Elattening  Method 
Kvith  the  similar  equation  obtained  from  the  Ray-Straightening  Method  shows 
that  the  two  equations  are  essentially  identical.  Moreover,  the  results  for 
the  transit  time  between  the  terminals  and  the  intensity  of  the  received  signals 
are  also  shown  to  be  essentially  the  same.  Thus,  the  predictions  for  the 
multipath  phenomenon  by  the  earth-flattenfcg- and  the  rav-straightening  methods 
were  shown  to  produce  the  equivalent  results  that  they  should  produce. 
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Two-Layer  Model  of  Line-of-Sight 
Refractive  Multipath  Propagation 


1.  INTRODUCTION 

The  theoretical  investigation  of  line-of-sight  (LOS)  refractive  multipath  prop¬ 
agation  has  not  attracted  much  attention  in  the  past  in  spite  of  the  experimental 
evidences  that  support  its  existence.  Perhaps  the  first  noteworthy  effort  on  this 
subject  can  be  found  in  the  work  of  Ruthroff.  1  He  assumed  that  the  earth's  atmos¬ 
phere  consists  of  several  concentric  layers,  each  of  which  has  a  constant  but  dif¬ 
ferent  gradient  of  the  index  of  refraction.  In  particular,  he  examined  a  two-layer 
atmosphere  in  some  detail.  This  model  has  the  obvious  advantages  of  physical 
and  mathematical  simplicity.  Here,  the  lower  layer,  which  is  immediately  above 
the  earth's  surface,  consists  of  the  standard  atmosphere,  while  the  upper  layer  is 
a  ducting  layer  of  sufficiently  large  thickness.  This  model,  although  overly 
simplified,  seems  to  deserve  a  more  detailed  study  because  it  helps  give  us  in¬ 
sight  into  the  physics  of  LOS  multipath  propagation.  Most  of  this  report  will  be 
devoted  to  drawing  inferences  from  this  model  in  a  systematic  fashion. 


(Received  for  publication  11  January  1983) 

1.  Ruthroff,  C.  L.  (1971)  Multiple-path  fading  on  line-of-sight  microwave  radio 
systems  as  a  function  of  path  length  and  frequency.  Bell  System  Technical 
Journal  50(No. 7):2375-2398. 


Several  years  after  the  appearance  of  Ruthroff's  work,  Pickering  and  DeRosa- 
placed  his  ideas  in  a  more  precise  mathematical  framework.  Instead  of  working 
in  the  actual  earth  coordinate  system,  these  authors  chose  to  transform  the  actual 
geometry  into  the  frame  of  reference  in  which  the  earth's  surface  is  flat.  This 
procedure  will  be  referred  to  as  the  earth-flattening  method  (EFM).  However, 

9 

the  work  of  Pickering  and  DeRosa"  is  incomplete  because  it  does  not  cover  all  the 
conceivable  relative  orientations  of  the  transmitter  and  receiver.  Moreover,  the 
two  orientations  they  considered  explicitly  had  to  be  treated  individually  because 
no  unified  set  of  equations  could  be  found  that  can  cover  both  of  them. 

A.  Malaga  and  S.  Pari2 3  of  Signatron,  Inc.  treated  the  problem  by  transform¬ 
ing  the  actual  situation  to  the  frame  of  reference  in  which  the  ray  trajectory  is 
piecewise  straight.  This  procedure  will  be  called  the  Ray -Straightening  Method 
(RSM). 

Since  the  EFM  and  RSM  represent  two  avenues  through  which  one  and  the  same 
phenomenon  can  be  examined  it  is  clear  in  principle  that  they  produce  exactly  the 
same  result  provided  that  they  are  followed  through  rigorously.  In  practice,  how¬ 
ever,  approximations  are  made  in  both  methods,  so  the  two  methods  cannot  be 
expected  to  lead  to  exactly  same  results,  although  they  should  be  almost  equiva¬ 
lent. 

In  this  report,  we  will  adopt  the  EFM,  remedy  the  incompleteness  in  Ref.  2 
mentioned  earlier,  and  then  show  the  equivalence  of  our  results  to  those  of  Ref.  3 
where  the  RSM  has  been  employed.  We  believe  that  establishing  such  an  equiva¬ 
lence  once  and  for  all  in  concrete  terms  will  be  a  worthwhile  endeavor.  It  should 
be  noted  that  one  can  observe  certain  connections  between  the  EFM  and  RSM  pre¬ 
dictions  even  before  doing  any  detailed  examination.  Specifically,  we  refer  to  the 
relation  that  geometrically  the  roles  of  the  ray  trajectory  and  the  layer  boundary 
are  interchanged  in  the  two  methods.  This  arises  from  the  fact  that  in  the  real 
geometry  the  ray  trajectory  segment  within  each  of  the  two  layers  may  be  approx¬ 
imated  by  an  appropriate  circular  arc  because  the  elevation  angle  of  the  ray  path 
remains  sufficiently  small  throughout  the  entirety  of  the  ray  trajectory.  Hie 
mirror  image  relationship  mentioned  above  leads  to  the  prediction  that  the  angles 


2.  Pickering,  L.  W.  ,  and  DeRosa,  J.  K.  (1979)  Refractive  multipath  model  for 

line-of-sight  microwave  relay  links,  1EEF  Trans,  on  torn  in .  Com -2  7  (No  -  3) 
1174-1182.  ~ . 

3.  Malaga,  A.,  and  Pari,  S.  (1981)  Theoretical  Analysis  of  Microwave  Propnga; 

tion.  Interim  Scientific  Report,  Contract  ^19628-80-T-oi<)<>,  Signatron,  ln<  . 
fTHe  findings  on  multipath  propagation  contained  in  this  report  are  included 
in  a  paper  by  S.  Pari,  "Characterization  of  multipath  parameters  for 
tropospheric  microwave  application"  submitted  to  1EE1.  Trans.  Antennas 
Propag. ) 
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of  take-off  and  arrival  in  the  two  methods  should  be  equal,  although  the  essential 
equality  of  the  delay  times  may  not  be  deduced  in  such  a  self-evident  manner. 

1  LINK  OK  SIGHT  REFRACTIVE  MULTIPATH  PROPAGATION 

We  assume  a  two-layer  system  in  which  the  lower  layer  is  the  standard  atmos¬ 
phere  while  the  upper  one  is  a  ducting  layer  with  a  refractive -index  gradient  less 
than  -l/a  *  -157  N  unit/km,  where  a  is  the  earth'  radius.  For  a  transmitter  (T  ) 
and  a  receiver  (Rx>  in  either  one  of  the  two  layers,  we  will  be  concerned  with  the 
study  by  the  EFM  of  the  propagation  of  signals  between  them.  As  was  mentioned 
in  the  Introduction,  the  work  in  Ref.  2  is  incomplete  in  the  sense  that  it  does  not 
cover  all  the  possible  orientations  of  the  terminals  relative  to  the  interface  between 
the  two  layers.  After  including  all  possible  relative  locations,  we  will  obtain  a 
unified  set  of  equations  to  determine  the  geometry  (and  hence  the  angles  of  take-off 
and  arrival)  of  a  ray  trajectory,  and  will  indicate  that  the  associated  trajectory 
angles  as  determined  by  the  EFM  and  RSM  are  essentially  identical. 

Then,  the  time  of  arrival  calculated  by  the  EFM  will  be  compared  with  that 
due  to  the  RSM  to  show  that  they  are  equivalent. 

Finally,  the  amplitude  factors  of  the  received  signal  calculated  by  the  two 
methods  will  be  shown  to  be  essentially  equivalent  also. 

2. 1  Angles  of  Take-Off  and  Arrival 

In  Ref.  2,  the  authors  assumed  that  the  height  of  Rx  is  not  smaller  than  that 
of  T  ,  and  considered  two  cases.  In  the  first,  T  is  in  the  lower  layer  and  R  is 

X  X  A 

in  the  upper  layer  (Condition  "A")  and  in  the  second,  both  terminals  are  found  in 
the  lower  layer  (Condition  "B").  Obviously,  this  is  not  a  complete  description  of 
the  problem.  Another  terminal  configuration  in  which  both  are  in  the  upper  layer 
(Condition  "C"),  must  be  added  to  the  system.  Completeness  is  also  contingent  on 
the  understanding  that  the  case  in  which  T  has  greater  height  than  R  can  be 
handled  by  applying  the  reciprocity  theorem.  In  this  report.  Conditions  "A  ', 

"B",  and  so  forth  in  Ref.  2,  will  be  designated  as  case  (A),  case  (B),  and  so  forth. 

2.  1.  1  UNIFIED  SET  OF  EQUATIONS  FOR  THE  NONDIRECT  PATHS 

Consider  a  transmitter  at  T  and  a  receiver  at  R  separated  by  a  distance  L 

X  A 

along  the  earth's  surface.  Let  the  heights  of  T^,  Rx>  and  the  layer  interface,  as 
measured  from  some  reference  level,  be  hT,  hR,  and  h0,  respectively,  and  let 


We  note  In  the  cases  of  practical  interest  that  the  elevation  angle  remains  suf¬ 
ficiently  small  throughout  the  ray  trajectory.  Under  such  a  circumstance,  the  ele¬ 
vation  angles  at  Tx  (the  take-off  angle)  and  at  Rx  (the  angle  of  arrival)  in  the 
EF  system  may  be  approximated  by  the  corresponding  angles  0^,  and  Sp,  respec¬ 
tively,  in  the  actual  geometry  (see  Section  A2.2).  Similarly,  the  elevation  angle 
®D  (the  grazing  angle)  in  the  EF  system  where  the  ray  hits  the  layer  interface  for 
the  first  time  after  it  leaves  Tx  may  be  approximated  by  the  angle  0q  in  the  actual 
geometry. 

Further,  let  the  lower  and  upper  atmospheric  layers  be  called  Layers  I  and  II, 
respectively. 

The  problem  to  be  solved  is  to  determine,  for  each  of  the  cases  (A),  (B),  and 
(C),  all  the  possible  ray  paths  that  connect  Tx  to  Rx.  We  will  show  that  this  can 
be  reduced  to  the  problem  of  determining  the  unknown  angles  0^,  0R,  and  for 
a  given  set  of  link  parameters.  However,  the  three  equations  for  them  will  be 
found  to  be  nonlinear  so  that  more  than  one  solution  for  the  problem  can  be  ob¬ 
tained.  The  collection  of  the  acceptable  solutions  that  are  picked  up  from  among 
these  general  solutions  form  what  we  will  call  the  multipath. 

However,  before  proceeding  to  write  down  the  unified  set  of  equations  for 
cases  (A),  (B),  and  (C),  a  few  remarks  are  in  order. 

First,  we  establish  a  sign  convention  for  the  angles.  In  the  following  treat¬ 
ment  of  the  problem  using  the  EFM  we  define  the  x  axis  (a  reference  axis)  to  be 
the  line  in  the  plane  of  the  ray  that  runs  parallel  to  the  earth's  surface  and  points 
positively  along  the  direction  from  Tx  toward  Rx-  Then  we  define  the  elevation 
angle  at  an  arbitrary  point  on  the  trajectory  as  the  angle  made  by  the  positive  tan¬ 
gent  at  that  point  according  to  the  standard  convention.  That  Is,  if  the  reference 
axis  must  be  rotated  counterclockwise  to  align  it  with  the  positive  tangent,  the 
elevation  angle  is  positive,  and  if  the  reference  axis  must  be  rotated  clockwise  to 
align  it  with  the  positive  tangent,  the  elevation  angle  is  negative.  When  the  eleva¬ 
tion  angle  at  a  point  on  a  circular  arc  is  transferred  to  a  central  angle  formed  by 
the  normal  to  the  arc  at  that  point  and  the  vertical  radial  of  the  circle,  the  sign 
and  magnitude  of  the  elevation  angle  reappear  conserved  as  the  central  angle  if  we 
regard  the  vertical  radial  as  a  new  reference  axis  and  assign  a  magnitude  and  a 
sign  to  the  central  angle  as  before.  By  following  this  standard  method,  which  is 
uniformly  applicable,  we  can  alleviate  the  confusing  sign-assigning  rules  used 
in  Ref.  2. 

Second,  since  the  ray  trajectory  within  a  layer  is  a  circular  arc  We  should 
be  awase  of  the  fact  that  a  statement  like  "a  terminal  with  a  given  height  lies  on  a 
ray  trajectory"  will,  in  general,  correspond  to  two  different  situations  where  the 
elevation  angle  of  the  ray  at  the  terminal  is  positive  or  negative.  That  this 


observation  is  of  importance  will  become  clear  when  we  start  writing  down  the 
equations  that  determine  the  ray  trajectory. 

Finally,  we  note  that,  depending  upon  the  relative  sizes  of  L,  d^,,  dp,  and 

L.j,  it  will  become  possible  for  a  ray  that  entered  the  second  layer  to  emerge 
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from  it  and  reenter  the  original  layers,  and  repeat  this  meandering  n  times  where 

n  =  0,  l,  2,  . . .  .  In  Figures  1  and  3  the  meandering  of  the  ray  trajectory  across 

the  layer  interface  is  represented  by  a  pair  of  short  vertical  dotted  lines  at  a  point 

where  the  ray  crosses  the  interface.  Referring  to  Figure  1  where  only  the  case 

of  t»P  <  0  is  sketched  for  illustrative  purposes,  we  have  for  case  (B) 

dp  =  R.(cos  0p  -  cos  Op) 

dT  *  Rf(cos  6^  -  cos  t>D) 

L  =  Rf(sin  tip  -  sin  0-j.)  +  n  •  2(R^  +  R)  sin  Op  +  2R  sin  Op 
+  Rj  (sin  Op  +  sin  Op) 


=  (2n  +  2)(Rf  +  R)  sin  Op  +  R^  sin  Op  -  R^  sin  0^, 


(1) 


In  the  above,  Rf  =  |Rjl  and  R  =  Rjj  where  R-,  i  =  1,  11,  are  the  radii  of  curvature 
Rt  -  I0fi/(L.  -  l/a)  of  the  ray  trajectory  in  Layers  I  and  II.  Note,  that  Eq.  (1)  is 
.  alid  for  0-p  S  o  and  6rI°  so  that  four  combinations  of  0^,  and  Op  are  possible. 
(These  possibilities  have  also  been  mentioned  in  Appendix  B  of  Ref.  3. ) 

By  introducing  new  variables 

I  =  sin  eR  .  n  =  sin  0-j,  ,  and  ?  =  sin  Op  ,  (2) 

Eq.  (l)  becomes 


?  -  n  +  t’1  ?  =  aL  ,  (3) 


4.  The  symbol  n  used  here  should  not  be  confused  with  the  refractive  index  n(r). 
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Figure  1.  Ray  Path  for  Case  (B) 


where 

°1 

=  dR/Rf 

"2 

=  dT/Rf 

°L 

II 

,r 

50 

°R 

=  R/Rf 

T-l 

=  (2n+2)(l  +  aR)  . 

Note  that  the  dimensionless  numbers  a  defined  in  Eq.  (4)  are  of  more  physical 
interest  than  those  numbers  introduced  in  Ref.  2.  Combining  Eq.  (3)  with  similar 
equations  for  cases  (A)  and  (C),  we  obtain  a  unified  set  of  equations  applicable  to 
all  cases, 

Vi  -  =  aj 

Vi  -  n2  -  Vi  -  ?2  =  ^ 


b|  +  crj  +  d?  =  a.  , 


(5) 


where  the  values  of  the  constants  a^.  ....  d  for  each  case  are  given  in  Table  1. 

(We  also  included  in  Table  1  the  values  of  constant  d  which  is  relevant  to  the 

o 

direct  paths  (see  Section  2.  1.  3). ) 


Table  1.  Values  of  Constants  for  Cases  (A),  (B),  and  (C) 


Case 

Co  ns  tant^'''-'\^ 

(A) 

(B) 

(C) 

al 

■°l/°R 

al 

a2 

°2 

a2 

-V°R 

b 

1 

c 

-1 

-1 

d 

a  *  =  (2n  +  1 )( 1  +  Op) 

t_1  =  (2n+  2)(1  +oR) 

-T-1 

d 

o 

CTo1  =  1  +°R 

0 

0 

Having  obtained  the  unified  set,  we  must  then  solve  it.  Since,  however,  the 
solution  procedure  for  the  unknowns  is  nearly  the  same,  we  will  discuss  only  the 
solution  for  n  as  an  example. 

It  will  be  seen  that  the  equation  for  n  that  can  be  derived  from  Eq.  (5)  is  of 
eighth  degree.  Unfortunately,  the  theory  of  equations  is  of  little  use  in  assessing 
the  nature  of  the  roots  of  such  an  equation.  On  the  other  fiand,  we  see  that  we  do 
not  have  to  solve  Eq.  (5)  exactly  because  the  practical  situation  we  will  be 
interested  in  is  the  one  for  which  both  of  the  conditions  r  a  a  and  cos  9  ~  1  are 
applicable. 

2.1.2  SMALL-  ANGLE  APPROXIMATION 

Since  we  have,  ir.  general,  a  series  expansion  for  cos  9  given  by 

cos  ®  =  1  -  j  sin“*  0  -  2~".  4  sin4  9  -  ••  • 


or  equivalently,  with  x  =  sin  9 


the  condition  cos  0  *  1  which  has  been  mentioned  above  indicates  that  the  square 
root  terms  in  Eq.  (5)  may  be  approximated  by  the  series  form  truncated  at  the 
second  term  (the  small -angle  approximation)  yielding 


b|  +  eq  +  d?  =  aL  .  (6) 

To  obtain  the  equation  for  r)  we  may  proceed  as  follows:  First,  from  the  first  two 
equations  we  get 

S2  -  r,2  =  -2(aj  -  ag)  ®  -2a.  .  (7) 

Next,  substituting  the  expression  for  S  obtainable  from  the  third  equation  into  the 
second  and  using  Eq.  (7)  we  obtain 

?  =  (Mg  rj2  +  M  j  p  +  MQ)/2b(cr7  -  a^),  (8) 

where 

M2  -  d2  -  (b2  +  c2) 

Mi  =  2oLc 

Mq  «=  2agd2  +  2a_b2  -a^  •  (9) 

Substituting  Eq.  (8)  into  Eq.  (7)  we  then  obtain  a  quartic  equation  given  by 
4 

(10) 

0 


where 


A.  =  M 
4  2  '  4b  c 


A,  *  4o,  c(M-  +  2d  ) 


(11) 


A2  -  2M2M0  -  4o2  (b2  -  c2)  +  8a_b2c2 
Aj  *  4aLc(MQ  -  4a_b2) 

A0  =  M5  +  8ft2La-b2  * 

It  is  to  be  noted  that  n  =  0  for  case  (A)  is  an  exception  where  we  have  A4  =  0  so 
that  the  resulting  equation  is  cubic  instead  of  quartic.  This  is  related  to  the 
apparent  discontinuity  of  Eq.  (l 0)  that  will  be  discussed  in  Section  2.  l.  4. 

We  think  that  this  represents  an  appropriate  stage  at  which  a  closer  compari¬ 
son  of  Eq.  (10)  with  the  similar  equation  derived  according  to  the  RSM  should  be 
made.  In  Ref.  3  Malaga  and  Pari  derive  a  quartic  equation  (their  Eq.  (2.42))  for 
an  unknown  x  =  that  is  applicable  to  all  conceivable  terminal  configurations. 
Rewritten  as  an  equation  for  8^.  it  is  given  by 

4 

J^A^  *  0  (12) 

0 

where 

\ 

A,  ■  (2a jb^  -  bjW 
A2  ■  (bl+Za^-o,)*2 

A1  '  Cblcl  •d2)*3 

A0  =  <cl  "  e2^4  •  (13> 

In  Eq.  (13).  $  a  D/Rq  corresponds  to  our  <»L  ■  L/R^  and  the  eight  constants  »j,  bj, 
Cj,  b2*  c2*  d2*  and  e2  dcPend  on  other  constants  ki2,  fcj.,  k^,  eT,  eR, 

R^t,  and  ReR  that  are  defined  in  terms  of  the  link  parameters.  The  intended 
comparison,  therefore,  requires  a  more  detailed  knowledge  of  these  constants. 
However,  considering  the  fact  that  Ref.  3  has  never  been  circulated  (it  was  intended 
only  for  RAOC/EE  internal  information)  it  makes  little  sense  to  go  into  detailed 
discussions  of  the  comparison  of  Eqs.  (10)  and  (13).  Therefore,  we  will  state 
here  simply  that  the  result  of  a  detailed  comparison  between  the  coefficients  shows 
that 


Ai6 


i  -  0.  1. 


4 
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where  6  *  a/Rj.  In  this  way  we  may  conclude  that  the  roots  for  n  ■  sin  $T  obtained 
by  the  EFM  are  identical  to  0-p  obtained  by  the  RSM.  This  shows  that  the  take-off 
angles,  the  acceptable  solutions  among  these  roots,  found  by  the  two  methods  are 
equivalent  to  the  extent  that  one  can  regard  sin  ■  ®T.  A  similar  statement  can 
be  made  concerning  die  angle -of -arrival.  It  should  be  noted  that  included  in 
Eq.  (3. 12)  of  Ref.  3  is  what  they  call  the  LOS  path  (which  is  the  same  as  the  direct 
path)  corresponding  to  the  cases  where  both  terminals  are  located  in  the  same 
layer  of  the  atmosphere.  In  contrast  to  this,  our  Eq.  (10)  is  applicable  to  the 
nondlrect  paths  only,  that  is,  to  the  cases  where  the  ray  crosses  the  layer  inter¬ 
face  at  least  once  and  hence  excludes  the  LOS  path.  Therefore,  strictly  speaking, 
equivalence  of  Eqs.  (10)  and  (12)  extends  only  up  to  our  nondlrect  paths.  However, 
the  equivalence  of  the  equations  for  the  LOS  path  in  Ref.  3  and  the  direct  path  in  our 
work  (see  Section  2. 1.  3)  can  be  readily  established  so  that  it  can  still  be  con¬ 
cluded  that  the  ray-path  equations  according  to  the  EFM  are  entirely  equivalent 
to  those  according  to  die  RSM. 

Incidentally,  a  quartic  equation  is  the  highest  degree  equation  that  can  be  solved 
algebraically.  The  sign  of  the  discriminant  indicates  the  kinds  of  roots  the  equa¬ 
tion  has.  Since  the  discriminant  is  a  function  of  the  link  parameters  via  die  co¬ 
efficients  of  the  equation  we  can  get  some  information  on  the  relation  of  the 
parameters  to  the  properties  of  the  roots.  In  practice,  however,  these  param¬ 
eters  enter  these  relations  in  such  an  intertwined  manner  that  no  useful  informa¬ 
tion  can  readily  be  extracted  from  these  relations. 

We  note  that  for  a  given  set  of  link  parameters  we  can  find  all  the  acceptable 
paths  that  interfere  with  die  layer  interface  in  the  following  manner:  We  consider 
the  solutions  to  Eq.  (6)  corresponding  to  nonnegative  integer  values  of  n  starting 
with  n  *  0  and  extending  up  to  a  certain  maximum  which  can  be  seen  to  exist  by 
examining  Eq.  (6).  For  each  value  of  n  the  four  roots  of  Eq.  (6)  must  be  exam¬ 
ined,  and  solutions  accepted  only  if  they  are  both  real  (that  is,  the  complex¬ 
valued  roots  must  be  discarded),  and  also  small  (that  is,  very  close  to  sero).  The 
collection  of  the  acceptable  solutions  of  Eq.  (6)  thus  found  will  constitute  the 
totality  of  the  paths  that  cross  die  layer  interface. 

2.1.3  DIRECT  PATH 

Although  it  was  not  made  unmistakably  clear  when  Eq.  (5)  was  first  written, 
it  was  implicit  that  it  was  intended  for  dealing  with  only  those  ray  paths  that  cross 
the  layer  Interface  (interface-crossing  (IFC)  paths). 

For  case  (B)  or  (C),  when  die  natural  propagation  of  the  waves  is  not  ob¬ 
structed,  there  can  exist  a  ray  path  that  is  not  included  in  Eq.  (5).  Such  a  path 
that  does  not  cross  the  layer  interface  is  usually  referred  to  as  the  direct  path. 


For  case  (A)  there  does  not  seem  to  be  any  obvious  criterion  according  to 
which  the  direct  path  can  be  defined,  except  for  the  requirement  that  it  should  be 
a  path  that  crosses  the  layer  interface  just  once;  that  is.  n  *  0  for  that  case. 
However,  it  is  possible  to  define  the  direct  path  for  case  (A)  in  the  following  way. 
Namely,  of  the  possible  three  solutions  to  Eq.  (6)  for  n  ■  0  of  case  (A)  one  can 
pick  up  that  real  solution  that  will  go  over,  for  the  fixed  link  parameters,  to  the 
direct  path  for  case  (B)  in  the  limit  when  the  height  of  the  layer  interface  is 
raised  to  approach  to  the  height  of  R^.  We  will  define  the  path  that  corresponds 
to  this  solution  to  the  cubic  equation  as  the  direct  path  in  conjunction  with  (or 
relative  to)  case  (B).  (Actually,  this  path  should  more  properly  be  regarded  as 
the  extension  to  case  (A)  of  the  direct  path  for  case  (B).  J  Similarly,  the  solution 
of  the  cubic  equation  that  goes  over  to  the  direct  path  for  case  (C)  in  the  limit 
when  the  height  of  the  layer  interface  approaches  to  that  of  will  be  defined  as 
the  direct  path  for  case  (A)  in  conjunction  with  case  (C).  A  fact  to  be  remembered 
in  this  definition  is  that  the  direct  path  for  case  (A)  is  a  solution  for  n  ■  0  of 
Eq.  (6)  for  case  (A). 

As  a  result  of  this  definition,  we  can  have  a  unified  set  of  equations  for  the 
direct  path  (see  Eq.  (14)  below]  with  the  understanding,  which  will  often  be 
repeated  in  what  follows,  that  die  appropriate  solution  for  case  (A)  must  be 
singled  out  from  the  three  possible  solutions.  Also,  it  should  be  possible  to 
establish  the  interconnection  among  our  nondirect  paths,  the  conventional  non- 
direct  paths  (where  the  direct  path  is  defined  only  for  cases  (B)  and  (C)] ,  and  the 
IFC  paths. 

W’ith  our  definition  of  the  direct  path,  Eq.  (14)  gives  the  exact  and  Eq.  (15) 
the  approximate  sets  of  equations  for  the  direct  path; 

Vl  -  c2  -  Vi  -  ?2 » 

Vl  -  n2  -  Vl  -  ?2  -  a2 

b€  +  cn  +  dQC  »  aL  ,  (14) 


€2  -  ?2  ■  -2a, 
n2  -  C2  •  -2a2 


bf  +  crj  +  d  C  ■  a 


(15) 


In  these  equations,  the  values  of  a^,  a^,  . . . .  and  dQ,  for  each  case,  are  listed 
in  Table  1. 

The  procedure  for  determining  the  direct  path  by  solving  Eq.  (IS)  Is  seen  to 
be  identical  to  that  for  the  IFC  paths  with  the  replacement  of  d  by  dQ. 

In  cases  (B)  and  (C),  however.  It  will  be  simpler  to  solve  Eq.  (15)  anew 
under  the  special  conditions  that  apply  to  these  cases;  that  Is,  b  ■  -c  and  dQ  «  0. 
This  leads  to  the  solution  given  by: 

n  /  2b  <*L  • 

With  b  *  1  and  -oR  for  cases  (B)  and  (C),  respectively,  we  may  combine  the  re¬ 
sult  as 


In  this  equation  we  set 


e*dT-dR 


-Rj  «  Rj  for  case  (B) 

+R  ■  Rtt  for  case  (C)  , 


where  Rj  and  Ry  are  the  radii  of  curvature  of  the  ray  trajectories  in  Layers  I  and 
II,  respectively. 

In  the  foregoing,  we  solved  Eqs.  (6)  and  (15),  the  small-angle  approximate 

versions,  for  both  die  nondlrect  and  direct  paths,  for  reasons  of  consistency. 

2 

In  contrast.  Pickering  and  DeRosa,  in  discussing  case  (B),  solved  the  exact 
relations  for  the  direct  path  (in  a  manner  equivalent  to  solving  our  Eq.  (14)] , 
while  solving  the  approximate  relations  for  the  nondlrect  paths.  Since  the  exact 
set  for  the  direct  path  for  cases  (B)  and  (C),  (Eq.  (14 )]  can  readily  be  solved 
analytically,  it  will  be  of  some  Interest  to  discuss  the  solution  of  Eq.  (14). 

For  cases  (B)  and  (C),  Eq.  (14)  leads  to  a  quadratic  equation  for  (  or  rj.  The 
solutions  are  given  by 


where  the  upper  (or  lower)  signs  off  and  n  go  together  to  be  consistent  with  the 
third  of  Eq.  (14). 

We  now  examine  the  problem  from  a  geometrical  rather  than  an  algebraic 
viewpoint  to  see  whether  or  not  both  possible  signs  represent  acceptable  results. 
Referring  to  Figure  2  for  case  (B)  for  example,  we  have 


sin  fa/  -  0-j.)  *  sin  (0R  -  <£>)  * 


where 


t  -  V  L2  +  e2  . 


Noting  that  sin  ■  e/t,  we  then  obtain 


0R  ■  sin”1  f  i  sin”1 


2^f  ’ 


From  this  we  get 


f 


2Rf  2Rf 


which  is  seen  to  coincide  with  the  result  in  Ref.  2  when  their  sign  convention  for 
is  taken  into  account.  The  result  given  by  Eq.  (19)  agrees  with  the  upper-sign 
solution  of  Eq.  (18)  for  case  (B).  Now,  a  question  remains  an  to  what  kind  of 
circumstances  would  give  rise  to  the  lower-sign  solution  of  Eq.  (18)  for  case  (B) 
and  whether  the  solution  is  physically  acceptable. 


Figure  2.  Direct  Path  (or  Case  (B) 


To  see  this,  suppose  that  we  are  given  a  problem  in  which  the  wave  can  travel 
along  the  circular  arc  of  radius  Rj.  If  we  consider  a  ray  which,  after  leaving  Tx, 
travels  along  the  dotted  circular  arc  of  radius  Rj  of  Figure  2  in  the  direction  indi¬ 
cated  by  the  arrow,  the  geometrical  treatment  similar  to  the  above  can  be  shown 
to  lead  to  the  solution: 


(20) 


This  is  Identical  to  the  lower-sign  solution  for  case  (B)  of  the  exact  set,  Eq.  (18). 
However,  we  can  readily  see  that  the  propagation  of  the  ray  from  to  along 
the  dotted  line  in  the  indicated  direction  is  clearly  in  violation  of  the  conditions  for 
the  ray  path  which  we  set  forth  in  the  beginning.  Therefore,  we  have  to  reject  the 
solution  given  by  Eq.  (20)  or  the  lower-signs  of  Eq.  (18),  and  will  be  left  with  only 
one  set  of  solutions  represented  by  Eq.  (19). 

Proceeding  analogously  we  can  show  that  the  lower-sign  solution  of  Eq.  (18) 
for  case  (C)  must  be  discarded  also. 

Combining  the  results  for  cases  (B)  and  (C)  we  see  that  the  exact  set  of  equa¬ 
tions  for  these  cases  yields  the  following  acceptable  solutions: 


—Id—  +  — £ — 

[M 

*•/  2ie  2\9L\  1 

i  \  *  / 

(21) 


We  should  note  that  normally  |2R|/t  1  so  that  the  second  term  on  the  right- 

hand  side  of  Eq.  (21)  is  approximately  equal  to  e/t  which  under  our  assumptions  of 


:  -  a  and  cos  H  =■  ;  may  bi  i  egarded  equal  to  e/L.  This  shows  that  the  exact 
result,  Eo.  (‘2l),  is  practically  the  same  as  the  approximate  one,  Eq.  (16). 

2.  1.-5  CONNECTION  FORMULAS 

According  to  the  definition  of  the  direct  path  given  in  the  previous  subsection. 
Eq.  (6)  with  n  *  0,  1,  2.  ...  will  provide  the  complete  description  of  path  geome¬ 
try  of  the  multipath  for  case  A;  that  is,  the  direct  plus  nondirect  paths.  However, 
for  case  (B)  or  (C)  it  will  not  become  complete  until  the  equation  for  the  direct 
path,  Eq.  (15),  is  brought  in.  Although  this  observation  is  sufficient  for  a  quali¬ 
tative  argument  it  is  not  sufficient  for  quantitative  description  because  the  con¬ 
tinuous  transition  of  the  multipath  phenomenon  between  cases  (A)  and  (B)  or  cases 
(A)  and  (C)  requires  a  little  more  precise  statement  which  we  might  call  the  con¬ 
nection  formulas  |see  Eqs.  (31)  and  (32)  below] .  (If  we  disregard  the  physical 
significance  attached  to  the  nonne^ative  integer  n  and  extend  its  range  to  include 
n  *  -1  on  purely  formal  grounds,  then  Eq.  (6)  corresponding  to  this  added  value  of 
n  is  seen  to  coincide  with  Eq.  (15)  for  cases  (B)  and  (C).  Then  Eq.  (6)  with  the 
extended  range  n  *  -1,  0,  1,  ...  will  describe  the  entirety  of  the  acceptable  paths 
for  cases  (B)  and  (C),  in  the  same  way  as  Eq.  (6)  with  n  *  0,  1,  2,  . . .  will  for 
case  (A).  ] 

To  determine  the  connection  formulas,  let  us  consider  the  transition  between 
cases  (A)  and  (B);  that  is,  the  situation  where  T^  is  located  in  Layer  I  while  Rx 
sits  exactly  on  the  layer  Interface.  Depending  upon  how  we  view  this,  namely, 
whether  we  regard  this  as  the  limiting  case  of  (A)  or  (B)  in  which  dR  —  0,  we  may 
have  to  employ  the  appropriate  version  of  Eq.  (6)  to  analyze  it. 

However,  before  dealing  with  the  above  problem  let  us  digress  a  little  to 
examine  the  special  case  of  (B),  called  (BS),  in  which  we  have  dT  *  dR  *  dg, 
because  it  will  help  us  gain  some  insight  into  the  handling  of  the  transition  case 
problem. 


Equation  (6)  for  case  (BS)  is  given  by 
f2  -  f2  =  -2os 
n2  -  f2  =  -2os 
?  -  n  +  t"1?  *  *L  , 

where  c*g  =  dg/R^.  This  may  be  reduced  to 


(22) 


(23) 
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where  the  upper  and  lower  lines  of  the  second  equation  correspond  to  ?  =  +n  and 
f  ■  -n.  respectively.  If  all  four  of  the  solutions  are  acceptable,  we  note  that  the 
two  solutions  with  C  *  -17  correspond  geometrically  to  those  ray  trajectories  that 
are  individually  symmetrical  with  respect  to  the  perpendicular  line  that  passes 
through  the  midpoint  between  the  terminals  (disregarding  temporarily  the  two 
short  dotted  vertical  lines  shown  in  Figure  3).  On  the  other  hand,  the  two  with 
C  *  +n  are  asymmetrical  individually  with  respect  to  the  perpendicular,  but  each 
is  the  symmetric  image  of  the  other  with  respect  to  the  perpendicular. 

By  calling  the  solutions  for  ?  =  +n  and  5  '  -n  classes  (a)  and  (b),  respectively, 
and  labeling  the  relevant  quantities  with  the  corresponding  subscript,  we  obtain 
the  following  quadratic  equation  for  each  class: 


t'2  „2  +  2og  r*2  -  a2  =  0 


(r‘2  -  4)  +  2«st‘2  -  a2L  =  0 


with  the  solution 


n  -  ■  ±  Ta . 


Figure  3.  Paths  for  Case  (BS) 


The  condition  for  the  roots  of  classes  (a)  and  (b)  to  be  real  are  given  by 
2  n  where 


"a  2(1  +  o„) 


Therefore,  the  necessary  condition  on  the  link  parameters  for  the  occurrence  of 
a  nondlrect  path  is  given  by  n  0.  Using  Eqs.  (26a)  and  (26b),  we  obtain 

3|  D 


(Rf  +  r  r 


for  classes  (a)  and  (b)  corresponding  to  the  upper  and  lower  lines,  respectively. 

When  nft  <  0  or  <  0  there  will  be  no  real  nondlrect  path  for  class  (a)  or  (b). 
On  the  other  hand,  when  n&  a  0  or  2  0,  additional  real  solutions  may  be  found 
corresponding  to  the  positive  integer  values  of  n  *  1,  2,  ...  up  to  In  1  or  N 
where  the  symbol  [q I  denotes  the  nonnegative  integer  closest  or  equal  to  but  not 
greater  than  q.  Since  it  can  be  shown  that  n^  -  nft  <  1  we  will  have  the  relation 
n^  <  n^  <  (n&  +  1)  which  shows  that  JnJ  can  exceed  by  unity  at  the  most. 

The  total  number  of  real  solutions  for  n&2  Ois  given  in  general  by  2  ( £nj  +  1). 
Using  this  information,  we  can  easily  determine  the  total  number  of  real  nondlrect 
paths  for  case  (BS),  which  consists  of  those  for  classes  (a)  and  (b). 

The  foregoing  shows  that  in  the  special  case  of  (BS)  the  quartic,  Eq.  (10)  can 
be  factored  into  two  quadratic  equations  thereby  enabling  us  to  find  the  analytic 
solutions  to  the  original  quartic  equation  in  a  very  simple  way. 

However,  as  soon  as  the  condition  dp  *  dR  is  relaxed  we  lose  the  symmetry 
relations  f  »  +n  and  ?  =  -rj  and  the  relevant  quartic  equation  can  no  longer  be 
factored  in  an  obvious  manner  as  was  possible  for  cue  (BS),  We  surmise  then 
that  we  will  have  to  resort  to  computers  to  solve  the  general  quartic  equation  for 
case  (B). 

We  can  now  profitably  return  to  our  original  question:  What  is  the  connection 
between  the  paths  for  the  limiting  case  of  (A),  called  case  (AI),  and  those  for  the 
limiting  case  of  (B),  called  case  (Bl),  where  Tx  is  found  at  the  same  position  in 
Layer  I  and  Rx  sits  on  the  layer  interface  ? 

For  cases  (AI)  and  (BI)  we  have  dR  =  0  so  that  Eq.  (6)  becomes 

?2  -  C2  -  0 
n2  -  52  =  ^ 

bf  +  cq  +  dg  *  .  (27) 

In  analogy  to  case  (BS)  let  us  continue  to  call  ?  *  +?  and  ?  *  -?  classes  (a)  and  (b), 
respectively.  Drawing  the  trajectories  for  each  of  the  classes  will  clarify  the 
geometrical  significance.  Proceeding  as  for  case  (BS),  by  referring  to  Table  1, 
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we  find  the  following  set  of  equations  applicable  to  classes  (a)  and  (b)  for  both 
cases  (AX)  and  <BI). 


2  -2  _  „ 
n  *  s  -  'Sag 

?  =  x  +  rj) 

where 

k  *(n)  *  d(n)  ±  b  = 
a 


Equation  (28)  leads  to  a  quadratic  equation 

(k  2  -  l)rj2  -  2 +  2o2k  ^  =  ®  (30) 


(28) 


o  *  +  »T 


T-1  ±  ! 


(2°}aR)s  i1 


(n)  for  case  (AI) 


+  (2n  +  2)  =  T^1(n)  for  case  (BI) 


with  solutions 


for  na  and  nb  for  cases  (AI)  and  (BI)  with  the  appropriate  expression  for  *  taken 
from  Eq.  (29). 

Here  again,  for  given  link  parameters  the  above  solutions  will  be  real  for 
such  nonnegative  integer  values  of  n  that  [nfl  b]  2  n  where 


The  necessary  condition  for  the  occurrence  of  the  nondirect  paths,  n&  b  2  0,  rna>' 
be  expressed  as  before  but  will  not  be  repeated  here. 

This  shows  that  the  quartic  equation  for  case  (AI)  or  (BI)  obtainable  from 
Eq.  (io)  is  expressible  as  the  product  of  two  quadratic  equations  (30)  for  classes 
(a)  and  (b)  of  each  case. 


From  the  foregoing  analyses  we  see  that  the  solutions  to  cases  (AI)  and  (BI) 
for  equal  n  are  not  exactly  the  same;  that  is,  not  continuous,  as  evidenced  by  the 
different  values  of  «  *  for  these  cases.  In  fact,  from  Eq.  (29)  we  see  that 

o^^n)  =  T^1<n)  for  n  =  0,  1,  2,  — 

a'^n)  =  t” *(n  -  l)  for  n  =  1,  2,  3 .  (31) 

a  a 

Here  we  note  for  class  (a)  that  the  partner  for  n  =  0  of  case  (AI),  a  *(0),  is  nnss- 

*  1  ^ 
ing  from  the  collection  of  for  case  (Bl). 

If  the  direct  path  has  been  counted  from  the  beginning  in  the  collection  of  the 
acceptable  solutions  for  case  (B),  then  in  the  limiting  situation  under  consideration 
we  will  have  the  direct  path  to  match  up  with  the  solution  corresponding  to  aa*(0) 
such  that  we  will  have 

for  n  =  0 

for  n  =  1,  2,  3, .  .  . 

(32) 

That  is.  the  cubic  equation  for  n  =  0  of  case  (AI)  which  factors  into  the  product  of 
a  linear  equation  corresponding  to  cr  *(0)  =  1  and  a  quadratic  equation  correspond- 
ing  to  e"  (0)  =  1  +  2qr  will  be  made  to  continue  to  tile  product  of  the  linear  equa¬ 
tion  for  the  direct  path  and  a  quadratic  equation  corresponding  toT^fO)  =  1  +  2»p 
for  case  (Bl).  Thus,  the  discontinuity  existing  between  the  solutions  for  cases  (A) 
and  (B)  of  Eq.  (6)  will  be  remedied  properly  to  achieve  continuity  by  considering 
for  case  (B)  not  only  the  solutions  to  Eq.  (6)  but  by  augmenting  it  by  direct 
path. 

A  discussion  analogous  to  the  foregoing  can  explain  the  situation  in  which  Rx 
is  situated  in  Layer  II  while  T^  is  found  on  the  layer  interface.  This  may  be 
regarded  as  a  limiting  situation  of  case  (A),  called  case  (All),  or  as  a  limiting 
situation  of  ease  (C),  called  case  (CII),  and  the  solution  in  the  multipath  phenom¬ 
enon  can  be  made  continuous  by  introducing  the  direct  path  for  case  (<■')  into  an 
otherwise  apparently  discontinuous  set  of  solutions  to  Eq.  (6)  for  cases  (AID  and 
(CII). 

In  concluding  the  subsection  let  us  comment  briefly  on  the  connection  between 
the  hidings  on  the  angles  as  given  in  lief.  3  and  in  this  report  (Sections  2.  1.  1- 
2.  i.  5).  As  we  have  seen  in  Section  2.  1.2  they  are  completely  equivalent  to  the 


Solution  corresponding 


The  limit  for  dp  -»  0  of 
direct  path  for  case  (B) 


Solution  corresponding 
to  r’^n  -  1) 
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extent  that  we  can  regard  sin  0  =  6  to  be  a  very  good  approximation.  However, 
the  manner  in  which  our  and  their  findings  are  presented  differ  considerably. 
Therefore,  it  will  be  of  interest  to  clarify  the  connection  between  the  two  results. 

In  what  follows,  we  will  attempt,  to  the  extent  possible,  to  do  that,  with  the  full 
understanding  that  Ref.  3  is  not  available  to  the  reader. 

Figures  2-15  and  2-18  of  Ref.  3  are  reproduced  with  some  modifications  in 
Figure  4.  These  modifications  include  the  addition  of  an  extra  dotted  or  solid  line 
to  the  single-lined  original  curve  and  attaching  arrows,  various  markings,  and 
numerals  for  the  integers  m  and  n,  to  facilitate  identification.  Here,  the  relation 
between  their  m  and  our  n  is  m  =  2n  +  1  for  case  (A)  and  m  =  2n  +  2  for  cases  (B) 
and  (C).  The  curves  in  the  figure  are  related  to  our  analysis  in  the  following  way: 

Consider  first  Figure  4(b)  where  the  curves  in  the  right  and  left  half -plane 
correspond  to  the  events  in  our  cases  (BS)  and  (CS),  respectively.  (Case  (CS)  is 
the  special  configuration  of  case  (C)  with  equal  terminal  heights;  that  is,  it  is  the 
counterpart  of  case  (BS).  ]  Confining  our  attention  to  its  right  half  for  the  moment 
we  can  see  that  the  curves  there  may  be  classified  into  two  distinct  groups.  The 
first  group  consists  of  a  family  of  parabolas  that  are  symmetric  with  respect  to 
the  abscissa  and  shrink  to  the  origin.  The  second  also  consists  of  a  family  of 
diminishing  parabolas  each  of  which  is  symmetric  with  respect  to  a  line  parallel 
to  but  a  certain  distance  from  the  abscissa.  Since  the  abscissa  of  the  figure  labels 
our  dg  =  RjOg  it  is  seen  that  the  parabolas  of  the  first  and  second  groups  represent 
the  angle -of -arrival  curves  corresponding  to  our  Eqs.  (24a)  and  (24b),  respectively, 
in  the  order  of  increasing  value  of  n.  At  the  height  of  the  layer  interface  indicated 
by  the  dashed  vertical  line,  the  values  of  the  angle  of  arrival  (according  to  the 
definition  in  Ref.  3)  for  n  =  0  for  classes  (a)  and  (b)  are  marked  by  the  symbols 
x  and  •  . 

In  Figure  4(a)  where  hR  -  h^,  *  0,  all  cases  of  (A),  (B),  and  (C)  occur.  For 
the  case  of  the  layer  at  the  height  of  Rx,  the  solutions  corresponding  to  classes  (a) 
and  (b)  of  Eq.  (30)  are  marked  by  x  and  •  as  before.  We  see  that  the  segment 
between  the  markings  x  and  ▲  of  the  doubly-lined,  outermost  curve  of  the  figure 
corresponds  to  n  =  0  in  our  notation.  The  various  portions  of  the  segment,  how¬ 
ever,  correspond  to  the  value  of  1  or  2  in  terms  of  the  integer  m  used  in  Ref.  3. 

The  inner  isolated  curve  (an  island)  of  the  figure  contains  a  segment  between  x 
and  A  that  behaves  in  the  same  way  as  in  the  outermost  curve  of  the  figure.  The 
segment  corresponds  to  n  =  1.  Now,  Ref.  3  gives  the  angle -of -arrival  curves  for 
hR  ”  hT  *e8S  than  m  chosen  *or  Figure  4(a).  Inspection  of  these  curves  shows 
that  the  number  of  isolated  islands,  each  surrounded  by  an  outer  island  in  suc¬ 
cession,  increases  as  h^  -  h.p  decreases.  In  the  series  of  diminishing  islands 
that  correspond  to  the  values  of  n  =  1,  2,  3,  ....  the  occurrence  of  the  same 
general  behavior  of  the  curves  that  was  observed  for  the  outermost  curve  (n  =  0) 
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is  apparent.  In  the  limit  of  the  value  of  hR  -  ttj,  vanishing,  which  is  the  case  for 
Figure  4(b),  the  islands  that  previously  appeared  isolated  now  have  points  of  con¬ 
tact  at  various  points  on  the  abscissa  of  the  graph.  A  comparison  of  Figures  4(a) 
and  4(b)  will  show  clearly  how  such  contacts  occur. 

2.2  Time  of  Arrival 

In  this  section  we  will  present  in  some  detail  the  result  for  the  time  of  arrival 
as  derived  from  the  EFM.  We  will  then,  as  in  the  previous  section,  compare  the 
EFM  result  with  the  result  from  the  RSM. 

2.2.1  NONDIRECT  PATH 

Figure  1  for  case  (B),  and  similar  diagrams  for  cases  (A)  and  (C),  show  that 
the  height  h  above  the  reference  level  of  an  arbitrary  point  on  the  ray  trajectory, 
where  the  corresponding  central  angle  (see  Section  2. 1. 1)  is  9(h),  is  given  by 

h  =  hD  +#£  (cos  0  -  cos  Sjj)  (33) 

where  hD  is  the  height  of  the  layer  interface  and 


Here  and  in  what  follows  he  (I),  for  example,  means  somewhat  loosely  that  the  point 
with  height  h  is  in  Layer  I. 

Now,  we  recall  that  the  modified  refractive  index  nm(r)  was  designated  by 

n  (h)  as  a  function  of  h.  Renaming  it  m(h)  we  sfee  that  it  may  be  given  by 
m 

m(h)  *  mjj  +  i^(h  -  h^)  (35) 

where 


In  the  above,  gj,  i  *  I  and  II,  are  the  constant  gradients  of  the  modified  refractive 
index  for  Layers  I  and  II,  respectively,  and  mD  »  m(D)  is  given  by 


mD  =  m(0)  +  gj  hD  . 


(36) 


Substituting  Eqs.  (33)  and  (36)  into  (35)  we  obtain 
m(h)  =  mD  ^(cos  6(h)  -  cos  ®D) 
Remembering  that 


(see  the  expression  for  R^(P')  somewhat  after  Eq.  (A26)J  the  above  becomes 

m(h)  *  M  -  cos  6(h)  (37) 

where 


M  =  nip  +  cos  0j-j 


Note  that  M  takes  on  a  different  value  corresponding  to  a  different  multipath, 
determined  as  sketched  in  Section  2.  1,  for  each  of  cases  (A)-(C).  To  see  that  M 
can  be  represented  in  different  forms  let  us  set  h  to  be  h  j  where  T  can  be  either 
T  or  R.  If  the  elevation  angle  at  J"x  is  called  ,  Eq.  (33)  gives  rise  to 


cos  eD  =  ^ 


<h. 


hp)  +  cos  0  j 


(38) 


which  when  substituted  into  Eq.  (37)  leads  to 


M  =  nrij-j  +  (h  j  -  hp)  +  cos  9-j.  =  m  y  +  cos  6y 


that  is, 

M  r  m^  +  cos  =  m  j  +  cos  0p  *  m^  +  cos  6p  (39) 

regardless  of  T  or  R  being  found  in  Layer  I  or  II.  The  equality  of  the  quantity 
in  +  cos  0  at  T  ,  R^,  and  at  the  point  where  the  path  crosses  the  layer  interface 
as  expressed  by  Kq.  (39)  is  of  wider  applicability.  That  it  is  actually  equivalent, 
under  the  small  angle  approximation,  to  Sr.ell's  law  can  be  seen  as  follows:  In  the 
EK  system,  Snell's  law  is  given  by 
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m(h)  cos  8(h)  *  const 


along  the  ray  path.  Rewriting  the  left-hand  side  as 
(m(h)  +  cos  8(h)J  +  (m(h)  -  lj  cos  8(h)  -  m(h) 
Snell's  law  becomes 


m(h)  +  cos  8(h)  *  const 


along  the  ray  path  under  the  assumption  that  m(h)  *  1.  We  see  that  what  is  shown 
in  Eq.  (39)  happens  to  display  the  approximate  constancy  along  the  ray  trajectory 
of  m  +  cos  8  at  the  three  special  points  on  it. 

To  calculate  the  transit  time  for  the  wave  to  travel  from  to  R^,  it  will 
prove  instructive  to  express  the  terminal  separation  L  and  the  curved  length  S  of 
the  path  as  integrals. 

Thus,  consider  the  expression  for  L,  Eq.  (1),  for  case  (B)  as  an  example. 
We  may  write 


L  * 


/  Rf  +  „ 


-8r 


f  (-R)  + 

8, 


D 


/  Rf 


-8, 


' D  R 

+  J  (-R)  +  J  Rp  )  cos  8  d8 


0T 


-8r 


-/  rf 


cos  8  d8 


T  -R 
x  x 


where  the  integration  over  8  Is  carried  out  from  T^  to  Rx  along  the  ray  trajectory. 
In  so  doing  we  noted  that  8  increases  in  Layer  I  whereas  it  decreases  in  Layer  II. 
it  will  be  seen  that  the  last  form  for  L  is  valid  universally  for  cases  (A)-(C). 

By  comparing  the  expressions  with  another  one  for  L  (see  the  third  of  Eq.  (5)| 


L  «  Rf  (b  sin  8p  +  c  sin  8^,  +  d  sin  8D>  (40) 

we  can  see  that  the  values  of  the  coefficients  b  and  c  as  listed  in  Table  1  can  also 
be  identified  with 


whereas  the  values  of  the  coefficient  d  do  not  lend  themselves  to  a  simple  exp  res 
sion  like  the  ones  in  the  above. 

In  a  similar  manner  it  can  be  shown  that  the  geometrical  length  S  of  the  path 
between  T^  and  R^  may  be  given  by 


S  -  -  f  dfl  «  Rf(bfiR  +  c<*T  +  d »D)  . 

Tx~Rx 

Returning  to  the  calculation  of  the  transit  time  r  for  the  link  we  see  that  it  is 
given  by 


f  ■  -  J  ^dfl/Ic0/m(®)J 


T  — R_ 
x  x 


where  m(8)  ■  m(h)  is  given  by  Eq.  (37)  and  cQ  is  the  velocity  of  light  in  vacuo. 
Therefore,  we  have 

r  -  ±  [m/*  (-se)  d  0-f  (-*f)cos8d8  -  i(MS  -  L) 

C°  1  T  -R  T  -»R  ° 

L  x  x  xx 

■  -jp£b(M  -  sin  ®p)  +  c(M  flT  -  sin  8^.)  +  d(M  8D  -  sin  8D)  J  . 


At  this  point  we  would  like  to  compare  the  above  results  with  those  found  in 
Ref.  3.  There,  the  terminal  separation  D  and  the  transit  time  are  shown  to 
be  given  by 


D  *  m<Rel  '  Rc2)flB  ‘  ReR®R  '  ReT®T 

cTrS  -  m(Rel  -  Re2)  sin  9fi  -  ReR  cos  9R  tan  9R  -  ReT  cos  »B  ten  «T 

where  c,  the  velocity  of  propagation,  represents  the  value  at  the  height  of  the 
layer  interface.  In  the  above,  m  *  0,  1,  2,  ...  denotes  the  number  of  crossing 
of  the  layer  interface  by  the  ray;  RgT  and  ReB  are  the  equivalent  earth  radii 
corresponding  to  the  layer  of  atmosphere  (Rg  j  for  Layer  I  and  *or  Layer  II) 
in  which  and  Rx»  respectively,  may  be  found;  9^,,  9R,  and  9R  are  the  angles 
of  take-off,  grazing,  and  arrival,  respectively. 

Comparing  Eqs.  (40)  and  (41)  with  the  expressions  for  D  and  cTrs  given  in  the 
foregoing  we  notice  that  they  are  parallel  in  the  sense  that  the  repeated  appearance 
of  the  same  coefficients  b,  c,  and  d  in  the  two  equations  of  our  set  are  copied  by 
the  appearance  of  the  same  coefficients  m(Rel  -  R^).  -  ReR.  and  -  ReT  in  the 
two  equations  of  their  set.  To  be  noted  also  is  the  appearance  of  9  in  D  at  the 
places  where  sin  9  appears  in  L.  This  is,  of  course,  a  reflection  of  the  fact  that 
roles  of  the  ray  trajectory  and  tee  layer  interface  are  interchanged  in  tee  EF  and 
RS  systems  (see  Appendix  B). 

To  compare  tee  results  in  more  detail,  we  note  the  following:  First,  our  Rf 
and  R  correspond  to  their  Re  j  and  -Re2 ,  respectively.  Next,  our  9T  and  9R 
correspond  to  their  9^  and  -9R,  respectively,  while  our  9D  corresponds  to  their 
9fi  for  cases  (A)  and  (B)  and  to  -9R  for  case  (C).  Therefore,  our 
Rjd  sin  9B  *  Rf  (2n  +  l)(l  +  oR)  sin  0D  for  case  (A)  and  Rf  d  sin  9D  *  ±Rf  (2n  +  2) 

(1  +  Op)  sin  9D  for  cases  (B)  and  (C),  n  *  0,  1,2,  ...  correspond  to  their 
m(Rel  -  Rej)®3  for  m  *  1,  3,  5,  ...  and  m  *2,  4,  6,  ....  respectively.  All 
combined,  it  shows  that  our  L  given  by  Eq.  (40)  is  approximately  equal  to  their  D 
on  the  condition  that  sin  9  •  9. 

What  should  be  of  more  direct  interest  to  us  is  to  see  whether  our  transit 

time  t  and  their  vRg  are  equivalent.  Numerically  they  cannot  be  equal,  of  course, 

because  they  are  defined  relative  to  different  values  of  propagation  velocity  cQ  and 

c.  A  direct  comparison  will  become  possible  when  the  transit  times  are  computed 

relative  to  the  same  velocity  of  propagation.  Thus,  for  example,  our  r  should  be 

compared  with  the  quantity  trsWco^  because  they  both- refer  to  the  same  velocity 

c  •  In  other  words,  we  should  examine  whether  tee  distances  c  r  and  ctoc  are 
o  o  RS 

equal.  That  this  is  essentially  the  case  can  be  seen  by  re -expressing  tee  three 
factors  of  tee  same  form  M9  -  sin  9  that  appear  in  Eq.  (4 1)  and  comparing  the 
result  with  the  expression  for  cTRg  given  in  Ref.  3.  Now,  using  Eq.  (39)  we  have 
for  the  first  factor  M9R  -  sin  9R  »  (mD  +  cos  9D>  sin  9R  -  sin  9R  «  cos  9Q  tan  9R, 
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and  similarly  for  die  second  factor,  while  for  the  third  factor,  we  have 
■  (M  *  l)  sto  *D  «  sin  With  the  above  we  see  that  cqt  agrees  with  cTpg 

showing  that  the  time  of  arrival  results  by  the  EFM  and  RSM  are  essentially 
equivalent  on  the  level  of  approximation  we  are  interested  in.  In  deriving  the 
above  result  we  noted  that  mD  *  1  and  8  *  sin  8  m  tan  8  are  valid  for  r  •  a  and 
cos  8  •  i,  respectively. 

If  we  define  the  time  delay  as  the  transit  time  relative  to  that  in  vacuo, 
then  it  is  given  by 

r  -  -t-  =  J_  (MS  -  2L)  (42) 

co  co 

The  corresponding  quantity  -  D/c  is  defined  as  the  relative  delay  in  Ref.  3. 
That  t  -  L/cq  s  0  and  *pg  -  D/c  £  0  can  naturally  be  understood  in  terms  of  the 
nature  of  the  EFM  and  RSM. 

At  this  point  we  would  like  to  make  some  comments  on  the  relative  delay  dia¬ 
grams  in  Ref.  3,  especially  in  conjunction  with  the  corresponding  angle -of -arrival 
diagrams.  Figures  2-25  and  2-28  of  Ref.  3  that  give  the  relative  delay  for  the 
paths  whose  angles  of  arrival  are  shown  in  their  Figures  2-15  and  2-18,  and  have 
been  discussed  in  Section  2. 1.4,  as  our  Figures  4(a)  and  4(b),  are  reproduced  here 
as  Figures  5(a)  and  505). 

First,  let  us  consider  Figure  5(b).  We  see  that  the  diagram  consists  of  a 
stack  of  figure  ^ ,  one  piled  on  top  of  another,  it  can  be  seen  that  the  bottom 
section  of  the  stack,  shown  in  double  lines,  corresponds  to  the  double-lined  curve 
in  Figure  4(b).  The  stacking  is  due  to  the  repetition  of  the  basic  pattern  of  the 
figure  for  n  ■  0  in  the  figures  for  n  ■  1,  2,  3,  . . . ,  which  corresponds  to  the 
similar  repetition  observed  in  Figure  4. 

Consider  the  time  delay  for  paths  with  n  ■  0  fet  some  positive  layer  height  in 
Figure  5(b).  Figure  4(b)  shows  that  there  are  two  paths  (marked  by  x)  belonging 
to  class  (a)  and  two  paths  (marked  by  •)  belonging  to  class  (b)  of  case  (BS).  Be¬ 
cause  the  paths  of  class  (a)  are  symmetric  to  each  other,  (see  Figure  3),  their 
relative  delays  are  identical;  they  correspond  to  one  and  the  same  point  of  the 
figure.  (In  Figure  5(b),  however,  we  deliberately  represented  them  as 
separated.  ]  The  same  applies  to  the  curves  for  n  *  1,  2,  3,  ...  of  Figure  5(b). 

In  Figure  5(b)  a  situation  analogous  to  Figure  5(a)  takes  place.  The  difference 
between  the  two  is  that  die  nonzero  value  of  hp  -  hT  in  Figure  5(a)  results  in  the 
spacings  between  the  successive  tiers  in  the  stacking  structure.  In  the  limit  of  a 
vanishing  value  of  hp  -  hT  these  spacings  reduce  to  zero  to  produce  Figure  5(b). 
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2.2.2  DIRECT  PATH 

The  expression  for  the  transit  time  ^  for  the  direct  path  is  given  by 


f  «Vd  •  L> 


where  Md  and  are  the  quantities  M  and  S  corresponding  to  the  direct  path. 

The  time  delay  for  the  direct  path  analogous  to  Eq.  (42)  is  given  by 

th  -  -2D  (43) 

d  c  c  d  d 
o  o 

If  the  experimental  circumstances  warrant  it,  one  may  be  interested  in  defining 
the  time  of  arrival  of  a  nondirect  path  as  the  transit  time  relative  to  that  for  the 
direct  path.  Then  the  time  of  arrival  ?  is  given  by 

r  *  r  -  r  .  _L(MS  -  M  .S  .)  (44) 

d  c  d  d 

o 

2.3  Amplitude  Factor 

In  this  section  we  will  find  the  power  (or  intensity)  factor,  and  from  it  the 
amplitude  factor,  associated  with  a  ray  as  it  is  received  at  Rx.  The  concept  of 
intensity  emerges  as  a  result  of  the  assumption  that  the  energy  which  was  initially 
confined  to  a  narrow  bundle  of  rays  that  emanates  from  T^,  and  centered  around 
the  ray  in  question,  continues  to  be  confined  within  that  same  bundle  throughout 
the  propagation;  the  increase  and  decrease  in  the  intensity  of  a  ray  occur  corre¬ 
sponding  to  the  contraction  and  expansion,  respectively,  of  such  a  bundle. 

Consider  a  small  bundle  of  rayB  centered  around  a  ray  with  take-off  angle 
0-p  and  angular  extents  of  60T  and  64  in  the  vertical  and  azimuthal  planes, 
respectively.  If  the  cross-section  and  the  intensity  of  the  bundle  are  6A(l)  and 
1(1)  at  a  point  with  unit  horizontal  distance  from  Tx  and  the  corresponding  quan¬ 
tities  are  6A(L)  and  I(L)  at  Rx,  then  we  have  6A(l)I(i)  *  6A(L)I(L)  according  to 
the  above  assumption.  When  the  elevation  angle  remains  small  throughout  the 
ray  path  it  can  be  shown  that 

6A(l)  *  60ji  •  64  •  cos  ©T 


6A(L)  «  |6L  •  L6#  •  sin  0  | 


where  6L  is  the  horizontal  extent  of  the  bundle  at  R^.  Then  in  the  limit  of  vanish' 
ing  cross -section  of  the  bundle  we  obtain  the  power  factor  (P  given  by 
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The  amplitude  factor  a  is  then  defined  by 


..4p  . 


The  corresponding  factors  for  the  direct  path  are 
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where  the  subscript  d  refers  to  the  direct  path. 

If  we  take  the  ratio  of  G*  to  ^d,  we  obtain  a  factor  that  approximately  repre¬ 
sents  the  power  of  a  nondirect  ray  relative  to  that  for  the  direct  ray.  It  is  given 
by 


t?a  lim  ^ 
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on  the  understanding  that  the  intensities  I^Cl)  of  the  direct  path  and  1(1)  of  a  non¬ 
direct  path  at  a  unit  distance  from  Tx  are  very  nearly  equal  to  each  other. 

When  the  free-space  propagation,  instead  of  the  direct  path,  is  chosen  as  the 
reference,  the  power  factor  (PQ  corresponding  to  Eq.  (47)  is  given  by 

/D  1JD  6A(l)  .  . 

o  1  lim  yiT  =  lim  6AlL)  1  -p  *  72 


where  the  subscript  o  refers  to  the  free-space  and  R  stands  for  the  slant  distance 
from  T^  to  R  .  The  power  factor  (PQ  relative  to  the  free-space  intensity  is  given 


<p„ 
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where  again  we  assume  that  1(1)  *  IQ(1).  This  is  the  power  factor  employed  in 
Ref.  3. 

We  note  that  the  power  factors  obtained  by  the  EFM  and  RSM  should  be  equiv¬ 
alent  since  they  are  determined  in  terms  of  various  angles  of  a  ray  path,  as  shown 
in  the  foregoing,  that  are  found  to  be  essentially  identical  (see  Section  2).  How¬ 
ever,  there  is  one  area  to  which  Ref.  3  directed  a  special  attention  while  Ref.  2 
did  not  touch  upon  it  at  all;  the  intensity  near  the  caustics.  We  observe  that  the 
authors  of  Ref.  3  discarded  the  RSM  which  they  had  followed  consistently  up  to 
this  point  and  opted  to  employ  the  EFM. 

As  is  well  known,  the  ray  theory  that  is  adopted  in  this  report  as  well  as  in 
Refs.  2  and  3  is  valid  only  under  certain  conditions.  In  particular,  the  intensity 
near  the  caustics  due  to  the  ray  theory  is  in  error  in  the  sense  that  it  grows  as 
one  approaches  a  caustic  and  becomes  infinite  at  the  caustic.  The  standard  method 
of  correcting  for  the  errors  is  to  examine  the  asymptotic  solutions  near  the  caustic 
of  the  wave  equation  and  compare  them  with  the  ray-theory  solutions.  This  has 
been  carried  out  in  Ref.  3,  obtaining  the  necessary  corrections  to  the  ray-theory 
results  near  and  at  the  caustics.  Since  their  work  seems  to  be  sufficiently  detailed 
and  accurate,  we  will  not  attempt  to  add  anything  more  to  it.  We  might  remark, 
however,  that  it  would  have  been  awkward  and  perhaps  rather  difficult  to  achieve 
the  same  results  had  they  continued  to  operate  in  the  RS  system. 

2.  3.  1  SMALL-ANGLE  APPROXIMATION 

This  is  the  approximation  in  which  the  nondirect  and  the  direct  paths  are 
determined  by  Eqs.  (6)  and  (15),  respectively. 

Evaluating  (Pas  shown  in  Eq.  (45)  by  using  Eq.  (6)  we  find  that 


1 

brjt  +  cfC  +  d$rj 


(50) 


The  corresponding  expression  for  is  obtainable  from  this  expression  by  the 

substitution  d  -*  d_  and  €  —  f  ,  ■  sin  and  so  on. 

o  d  _  n  , /„ 

Therefore,  the  amplitude  factor  a  *  (P'7  w^)1'  can  be  calculated  as  soon  as 
€,  n,  and  5  for  the  nondirect  and  direct  paths  are  found  by  solving  Eqs.  (6)  and 
(15),  respectively. 
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For  the  direct  path  for  case*  (B)  and  (C)  [but  not  for  caae  (A))  where  apecial 
conditiona  dQ  »  0  and  b  *  -c  hold,  the  expression  within  die  abeolute -value  algn  of 
Eq.  (50)  reducea  to  a  *  R^/L  so  that  we  obtain  ■  L*2.  This  ahows  that  (P^ 
ia  approximately  equal  to  (PQ.  establishing  the  corresponding  equality  between  the 
relative  power  factor  G>  in  this  report  and  (PQ  employed  in  Ref.  3. 

2.3.2  EXACT  RESULT 

It  may  be  recalled  that  in  Ref.  2  the  amplitude  factor  a  was  calculated  for 
case  (B)  by  employing  the  approximate  solution  for  the  nondirect  path  and  the  exact 
one  for  the  direct  path.  This  apparent  nonuniformity  should  not  be  a  source  of 
any  special  concern  since  our  assumptions  that  r  *  a  and  cos  6  =  1  for  the  ray 
path  will  guarantee  that  the  exact  and  the  approximate  solutions  for  the  path  are 
practically  the  same  after  all.  [See  the  end  of  Section  2. 1.  3.  ] 

Nonetheless,  for  the  sake  of  completeness  we  will  give  here  the  expression 
for  (?  derivable  from  the  exact  set  of  Eq.  (5).  It  is  given  by 

_ g  V  1  -  n2 _ 

brjg  V  1  -  I2  +  cg|  >/l  "  I2  +  din  V  1  -  g2 

This  expression  for  (P  goes  over  to  that  in  Eq.  (50)  when  the  square -root  terms 
are  approximated  by  unity.  The  exact  expression  for  may  be  obtained  from  the 
above  by  the  method  mentioned  in  the  previous  section. 
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3.  CONCLUSION 

In  this  report  the  two-layer  model  of  the  line-of-sight  refractive  multipath 
propagation  for  microwave  links  has  been  examined  in  accordance  with  the  ray 
theory. 

To  describe  the  phenomenon  we  employed  the  earth-flattening  coordinate  sys¬ 
tem.  The  result  obtained  was  compared  with  the  corresponding  work  found  by  the 
workers  of  Signatron,  Inc. ,  who  adopted  the  ray-straightening  system.  Specifi¬ 
cally,  the  angles  of  take-off  and  arrival,  the  time  of  arrival,  and  the  amplitude 
factor  were  examined  in  detail. 

From  the  purely  theoretical  point  of  view,  the  results  for  these  quantities 
should  not  depend  on  which  of  the  two  methods  is  used.  However,  the  equivalence 
of  the  results  is  not  self-evident  until  it  is  shown  in  concrete  terms. 

It  was  shown  explicitly  in  this  report  that  the  relevant  angles,  the  time  of 
arrival,  and  the  amplitude  factor  determined  by  the  earth -flattening  and  ray¬ 
straightening  methods  are  in  fact  essentially  equivalent. 
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Review  of  ■  .ay  Theory 


Al.  RAY  THEORY  IN  GENERAL 

We  will  summarize  here  some  well  known  pertinent  results  of  ray  theory. 
Given  a  linear  homogeneous  partial  differential  equation 

f(V.  3tH  =  0  ,  (Al) 

where  V  is  the  gradient  operator  and  3^  =  3/3^,  we  get  a  primary  solution  of  the 
form 

1  =  g(^)  , 


=  a  •  x  +  a4t  (A2 ) 

with  vectors  a  =  (a^,  a2»  a3>  and  x  =  (Xj.Xg.Xg),  where  the  coefficients  ai  satisfy 
f(a,  a^)  =  0  . 

The  solution  given  above  also  satisfies  the  so-called  equation  of  characteristics  of 
the  original  differential  equation  defined  by 
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f(vt,atti)  =  o  . 


If  we  specialize  to  the  function 
f(b,b4>  =  b2  -  b2/c(x) 
we  obtain 


f(V,  9t>C 
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V~t 
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which  represents  the  scalar  wave  equation  when  c(x)  is  the  speed  of  propagation  at 
the  point  x.  The  primary  solution  to  this  equation  may  be  given  by  Eq.  (A2)  with 

0  =  a  •  x  -  c  t 
—  —  o 

where  a  =  d/ldl  with  the  direction  cosine  vector  d  =  (d-.d^dg),  and  cq  is  the 
speed  of  propagation  in  the  empty  space.  When  a  more  general  form  for  0 , 

0  =  S(x)  -  CQt  (A  5) 

is  considered,  the  equation  of  characteristics  becomes 

(VS)2  =  n2  (At>) 

where  n  =  cQ/c  is  the  refractive  index.  This  equation  is  commonly  referred  to  as 
the  eikonal  equation,  and  is  of  fundamental  importance,  leading  to  the  concept  of 
rays.  Because  of  the  choice  made  for  0  in  Eq.  (Aii)  it  will  be  seen  that  the  func¬ 
tion  t i  =  g(S(x)  -  c^t)  constructed  in  terms  of  the  exact  solution  to  the  eikonal 
equation  cannot  satisfy  the  original  differential  equation,  that  is.  the  wave  equation. 
Noting  first  that  the  solutions  in  terms  of  sinusoidal  motion  is  our  basic  interest 
and  second  that  the  amplitude  of  the  oscillation  should  be  introduced  as  an  unknown 
function  of  spatial  coordinates,  we  are  led  to  consider  the  solution  to  Kq.  (A4) 
given  in  the  form 

i;  =  A(x)  exp  ( -iu)0  (x,  t)./c^|  =  A(x)  exp  Inct  -  ikS(x)| 

where  k  -  u/c^.  Substituting  the  above  into  Icq.  (A 4)  and  separating  tin  real  an 
imaginary  parts  we  obtain 
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For  S  to  be  a  solution  of  the  eikonal  equation  the  last  term  of  the  first  equation 
must  be  zero  which  means  that  the  frequency  must  be  infinitely  large.  When  S  is 
a  solution  to  the  eikonal  equation  it  will  serve  as  a  good  approximation  to  the  solu¬ 
tion  to  the  wave  equation  provided  that  the  second  term  in  the  second  of  Eqs.  (A7) 
is  small  in  comparison  to  the  first.  Note,  however,  that  it  is  not  sufficient  to  say 
that  k  must  be  large.  We  have  to  specify  how  large  k  must  be  relative  to  the 
physical  conditions  of  the  problem  under  consideration.  The  condition  we  are 
looking  for  can  be  found  through  order  of  magnitude  considerations,  and  may  be 
stated  as  follows: 

A  solution  to  the  eikonal  equation  will  serve  as  a  good  approximation 
to  the  phase  function  in  the  solution  to  the  wave  equation  if  the  frac¬ 
tional  change  in  the  velocity  gradient  over  a  wavelength  is  small 
compared  with  the  velocity  gradient  itself. 

Returning  to  the  question  of  how  to  construct  the  solution  to  the  eikonal  equation, 
we  can  show  that  it  is  recognizable  as  a  problem  of  finding  a  function  with  a  pre¬ 
scribed  magnitude,  but  not  direction,  of  the  gradient.  It  has  been  known  that  the 
solution  leads  to  the  concept  of  the  rays  and  the  wavefront,  giving  rise  to  the 
vector  form  of  the  ray  equation  for  the  general  case  which  reduces  to  simpler  forms 
for  certain  special  classes  of  the  refractive  index.  Moreover,  it  is  known  that 
these  ray  equations  can  also  be  derived  from  Fermat's  principle. 
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where  dl  is  the  line  element  of  the  ray  trajectory,  as  follows: 

(1)  General  case:  n  =  n(x): 

Representing  the  ray  path  by  x  =  x(s)  with  a  parameter  s,  we  get  a  line 
element  given  by  dl  =  ^  x'(s)2  ds.  Then  Fermat's  principle  is  cast  into 

2 

6  y  I.lx(s),  x'(s)|  ds  =  () 
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with  the  Lagrangian  L  =  n(x(s)j  ^  x'(s)^.  The  Euler  equation  corresponding  to 
this  leads  to  the  vector  form  of  the  ray  equation  given  by 


CA8) 


(2)  The  gradient  of  n  is  always  in  the  radial  direction: 

This  corresponds  to  the  case  where  n  is  spherically  stratified: 
n  =  n(r),  r  =  |x| .  If  at  a  point  on  the  ray  trajectory  its  local  elevation  angle  is  9 
and  its  position  vector  makes  an  angle  <t>  with  a  reference  direction  through  the 
center,  we  have  dl  =  V  1  +  (r^')2  dr  and  L($(r),  $'(r)J  =  n(r)  V  1  +  (r$')2  so  that 
the  resulting  Euler  equation  reads 


rn  cos  9  =  const 


(A9) 


(3)  The  gradient  of  n  is  constant  in  direction: 


Describing  n  as  n(h)  in  terms  of  the  rectangular  coordinates  (h,  s)  we  get 
dl  «  V  1  +  s'2  dh  and  L|s(h),s'(h)|  *  n(h)  yj  1  +  s'2.  The  Euler  equation  is  given 

by 


n  cos  0  =  const  .  CAlO) 

This  relation  is  generally  called  Snell's  law.  Because  of  its  relation  to  Eq.  (A  10) 
the  relation  expressed  in  Eq.  (A8)  is  usually  referred  to  as  the  generalized  form 
of  Snell's  law.  We  might  mention  that  the  name  "Snell’s  law"  is  frequently  applied 
also  to  the  relation  given  by  Eq.  (A9),  although  some  authors  may  prefer  to  call 
it  the  Bouguer'  formula. 

A2.  APPLICATION  OF  RAY  THEORY  TO  A  SPHERICALLY 
STRATIFIED  ATMOSPHERE 

The  two-layer  model  of  LOS  refractive  multipath  that  we  want  to  investigate 
assumes  a  spherically  stratified  atmosphere  n  *  n(r)  in  a  spherical  coordinate 
system  having  its  origin  at  the  earth's  center.  Then  the  relation  given  by  Eq.  (A9) 
prevails  in  the  plane  that  includes  the  ray  trajectory.  In  this  section  we  will  first 
examine  certain  aspects  of  the  LOS  propagation  over  the  spherical  earth  for  a 
general  spherically  stratified  atmosphere.  Then  the  same  situation  will  be  dis¬ 
cussed  from  the  viewpoint  of  the  EFM  and  the  RSM.  These  results  will  finally  be 


3.  See,  for  example.  Born,  M.  and  Wolf,  E.  (1965)  Principles  of  Optics, 
Pergamon  Press,  New  York,  p.  123.  “  ~~~ 


specialized  to  the  case  of  practical  interest  where  the  gradient  of  the  refractive 
index  is  constant  and  the  transmitter  and  the  receiver  are  found  near  the  earth's 
surface.  Figure  Al  depicts  schematically  what  happens  in  each  of  three  situations 
mentioned  above.  Figure  Ai(a)  shows  the  propagation  of  waves  in  the  real  world, 
that  is,  over  a  spherical  earth,  while  Figures  Al(b)  and  Al(c)  sketch  the  same 
event  as  seen  in  the  EF  and  RS  frame  of  references. 

A2.1  Propagation  Over  a  Spherical  Earth 

Referring  to  Figure  Al(a)  consider  a  ray  which,  starting  from  a  point  P  with 
polar  coordinates  (r  1#  *-/2),  propagates  through  a  spherically  stratified  atmosphere 
n  *  n(r)  and  passes  through  another  point  P  with  coordinates  (r, ?).  When  the 
geometrical  relation 

_  dr  cot  6 
de  ■  — - -  . 

where  4  =  t/2  -  4  and  0  is  the  local  elevation  angle,  is  combined  with  the  ray  path 
equation  for  this  case  *  compare  with  Eq.  (A9)J 

r  n  cos  0  =  A 

we  obtain  the  following  differential  equation  for  4: 

4'(r)  ■  =  - -  -A  -  -  .  (A  11) 

r  V  (rn)2  -  A2 

From  this  equation  we  obtain  the  solution 

r 

♦(r)  =  A  f  - —  dr~  ~ 

r  V (rn)2  - 

as  well  as  the  relation 


*"<r> 


2(rn)2  -  A2  +  r3  nn 
2  [(rn)2  -  A2]3'2 
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Next,  let  us  define  the  radius  of  curvature  R(P)  at  a  point  P  on  the  ray  trajectory, 
as  expressed  by  y  *  y(x),  by 


IMP) .  . 

y"(x) 


with  a  minus  sign  in  front.  It  follows  from  this  definition  that  the  radius  of  curva¬ 
ture  so  defined  is  an  algebraic  quantity  such  that  it  will  be  <0  (or  >0)  for  a  tra¬ 
jectory  which  is  concave  up-  (or  down-)  ward.  In  terms  of  the  polar  coordinates 
(r,  $  ■  wl 2  -  ♦)  it  becomes 

(1  +  r2d'(r>2)3/2 

RIP)  ,  + - - 

2*’<r)  +  rf '(r)  +  r  ♦Mr)4 

which  on  substitution  of  Eqs.  (All)  and  (A12)  gives  rise  to6 


n'  cos  6/n 


(A13) 


A2.2  Ewth-FUtteninjr  Transformation 

The  trajectory  of  the  ray  in  the  earth-flattening  coordinate  system  can  be  found 
in  the  following  way.  Imagine  an  observer  who  moves  along  the  earth's  surface 
and  measures  the  distance  traveled  s  and  the  height  above  the  ground  h  of  the  ray 
trajectory.  Based  on  these  data,  construct  a  rectangular  coordinate  system  (s,h). 
Then  the  curve  which  describes  the  loci  of  h  as  a  function  of  a,  that  is,  h  *  h(a), 
represents  the  ray  trajectory  in  the  EF  system.  Analytically,  the  EF  transforma¬ 
tion  is  characterized  by  the  relations  as  given  below : 

s  *  ad 


h  *  r  -  a 


(A  14) 


The  transformed  ray  trajectory  in  the  EF  system  is  sketched  in  Figure  Al(b)  in 
which  points  A  and  P  of  the  actual  configuration  are  mapped  to  points  A  and  P', 


6.  This  relation  was  derived  geometrically  in  Millington,  G.  (1957)  The  concept 
of  the  equivalent  radius  of  the  earth  tropospheric  propagation,  Marconi  Rev. 
20(No.  l26):79-93. 


respectively,  where  the  point  P'  has  the  coordinates  (s.h).  The  radius  of  curva¬ 
ture  Rj(P')  at  P*  of  the  transformed  trajectory  is  given,  in  accordance  with  the 
sign  convention  adopted  earlier,  by 


+  h'(sr 


R,(P')  - - : 

1  h"(s) 


Noting  that 


dh/dr  j 
ds/dr  a$'(r) 


.  t  _  d  dh  _  ♦"(*“) 

h'*»  -  d?  31  -  -  ^^3 


and  using  Eqs.  (All)  and  (Al2)  we  obtain 


R  j(P')  =  + 


[l  +  aV(r)2] 
a*"(r> 


r2[l-{2(.-,2/r2>]: 

a€<2  -  e2  +  rn'/n) 


where 


C(r)  *  A/rn(r)  *  cos  0(r)  , 


At  this  point  we  want  to  make  a  few  remarks  regarding  Eq.  (A  16).  By  substituting 
n'  *  rn2/AR(P)  from  Eq.  (Ai3)  into  Eq.  (Al6)  we  may  express  l/R^P')  in  the  form 


1  1  _  1  1 
R , (P*>  "  ®  R(Phr„(r)  "  av"(r)  * 


where 


vp(r)  «  -  [l  -  C2(l  -  a2/r2)]3  2 

r2  [i  -  52(.  -  .V)]  m 

v  <r )  «  - ^ = - i - 


a2C(2  -  f2) 


In  the  case  of  practical  interest  where  r  *  a  and  cos  0*1  throughout  the  ray 

trajectory  we  further  note  the  following. 

First,  v0(r)  *  1  and  v  (r)  *  1  under  these  conditions  so  that  the  above  relation 
r  a 

among  the  curvatures  simplifies  to 

R^P')  "  i  *  R^Pl  "  a"  •  (Al7) 

The  relation  shows  that  the  relative  curvature  between  the  ray  trajectory  and  the 
earth's  surface  taken  at  the  properly  corresponding  points  will  be  conserved  in 
the  transition  from  the  actual  to  the  EF  geometry.  As  already  seen  above 
]Eq.  (A  16)],  in  the  general  case  where  the  stated  conditions  need  not  be  satisfied, 
the  relative  curvature  in  the  EF  system,  that  is,  the  difference  between  the  curva¬ 
tures  l/RjiP')  of  the  trajectory  at  P'  and  l/«  of  the  earth's  surface  at  A,  does  not 
equal  the  relative  curvature  in  the  actual  system,  that  is,  the  difference  between 
l/R(P)  of  the  ray  trajectory  at  P  and  l/a  of  the  earth's  surface  at  A.  It  should  be 
realized  that  the  validity  of  the  conservation  depends  on  the  fulfillment  of  the 
conditions  r  *  a  and  cos  0  *  1. 

Second,  from  Eq.  (Ai6)  it  is  seen  that  we  will  have  Rj(P')  ^  0  for  $"(r)  ^  0. 
Under  the  stated  conditions  this  latter  relation  corresponds,  from  Eq.  (12),  to 


-n'  5  i  *  157  N  unit/km  . 

Before  closing  the  subsection  we  will  make  brief  comments  on  the  concept  of  the 
modified  refractive  index  which  plays  a  significant  role  in  the  treatment  of  various 
problems  according  to  the  EFM,  and  also  on  the  connection  between  the  elevation 
angles  in  the  actual  and  the  EF  coordinate  systems.  The  exact  slope  at  P'  of  the 
trajectory  in  the  EF  system  is  given  by 


dh  1 
ds  *  a4  (r) 


J"L(h>2  -  A2 

j.  nr - 

X 


where 


n  (r)  *  -  Mr)  ■  (l  +  jM  n(h)  ■  n  (h) 
m  a  \  a  f  m 


(A  18) 


(Al9) 


defines  the  modified  refractive  index  with  n(r)  ■  n(h)  and  A  *  A/a.  Therefore,  in 
terms  of  the  elevation  angle  6(h)  at  P'  we  have 


tan  0(h)  * 


V  •  *2 

K 


If  our  attention  is  confined  to  the  case  where  h/a  «  1,  the  above  relation  becomes 
equivalent  to 

n  (h)  cos  0(h)  *  X  •  const  . 
m  ^ 


Recall  from  Eq.  (A  10)  that  this  represents  the  ray  trajectory  equation  for  a  med¬ 
ium  with  refractive  index  n  m(h)  as  given  by  Eq.  (A  19).  To  determine  the  relation 
between  the  elevation  angles  0(r)  in  the  actual  geometry  and  0(h)  in  the  EF  system 
we  note  that  the  slope  in  the  former  is  given  by 


tan  0(r) 


dr  _  a  dh 
r  d^  r  ds  * 


where  use  has  been  made  of  Eq.  (Al4).  Since  dh/ds  is  the  slope  of  the  trajectory 
in  the  EF  system  and  hence  is  equal  to  tan  0(h)  we  obtain  the  desired  connection: 

tan  0(h)  ■  £  tan  0(r)| 

a  Ir-a+h 

Under  the  conditions  of  practical  interest  we  may  make  the  approximation  that  the 
factor  r/a  is  unity  by  dropping  the  contribution  from  h/a.  leading  to  the  approxi¬ 
mation  0(h)  ■  0(a  +  h).  This  means,  in  particular,  that  the  take-off  and  arrival 
angles  determined  experimentally  may  be  identified  with  the  corresponding  angles 
that  can  be  calculated  in  the  EF  system.  We  should  note,  on  the  other  hand,  that 
the  factor  r/a  in  Eq.  (A  19)  must  not  be  approximated  by  unity.  This  is  because 
the  propagation  phenomena  could  depend  sensitively  on  a  relatively  small  change 
in  the  refractivlty  so  that  it  also  becomes  necessary  to  retain  the  contribution 
from  h/a  even  though  this  number  may  be  very  small  in  comparison  to  unity. 

A2.3  Ray-Straightening  Transformation 

The  transformation  from  the  actual  geometry  to  the  RS  system  can  be  carried 
out  in  a  manner  analogous  to  the  EF  transformation.  Namely,  we  imagine  an 
observer  who  measures  the  distance  travelled  I  and  the  depth  to  the  earth's  sur¬ 
face  p  by  moving  along  the  ray  path  in  the  actual  geometry.  When  a  rectangular 
coordinate  system  (i,p)  (shown  in  solid  curves  in  Figure  A  1(c)]  is  constructed  it 


6 


describes  the  event  in  the  RS  coordinate  system.  Consider  a  point  Q'fc.p)  on  the 
transformed  curve  for  the  earth's  surface  and  call  the  radius  of  curvature  there 
R2(Q').  The  radius  of  curvature  can  be  found  in  a  way  parallel  to  the  EF  case. 
First,  from  Figure  Al(a)  we  note  the  geometrical  relation 

di  *  dr/sin  0 
from  which  we  obtain 

l'(r)  *  -  - -  .  CA20) 

V( rn>*  -  A2 


f"(r) 


A2(rn)’ 


<rn>2  -  A2 


and 


f(r)  * 


/ 

ri 


dr  rn 

V  (rn)2  -  A2 


(A21) 


Next,  we  require  the  expression  for  p(r).  The  simplest  way  to  find  this  is  to  note 
that  Z  OPQ  =  6  from  which  one  can  obtain  directly  that 


p(r)  =  r  cos  0 


-  r2  sin2  0  . 


(A22) 


To  determine  R2 (Q')  let  us  consider  the  mirror  image  [shown  in  the  chained- 
lines  in  Figure  A  1(c)]  of  the  RS  picture  outlined  above.  If  the  same  letters  t  and  p 
are  used  for  this  new  system,  the  radius  of  curvature  R2"(Q")  at  Q"  will  be 
given,  in  accordance  with  our  sign  convention,  by 


R2"(Q")  -  - 


[.w]3/2 

p"(i) 


47 


Since  we  have 


p"(i)  -  l  "(rfrp^r) 

"  i'(r>3 


where  the  prime  refers  to  the  derivative  with  respect  to  the  argument,  and 
obviously  R2(Q')  =  -R2"(Q"),  we  obtain 


(i'2  +p'2)3/2  1  +  (1  -  C2)p*2] 


3/2 


R?(Q')  =  +  #.ii  ,-m  =  9  9 

2  ip  ~pi  (1  -  €2>p"  +  ?2p'(rn)'/rn 


CA23) 


On  the  right-hand  side  of  this  equation,  I  and  p  are  considered  as  functions  of  r, 
and  use  has  been  made  of  Eqs.  (A20)-(A22).  As  before,  using  Eq.  (Al3)  we  obtain 
the  following  relation  from  Eq.  (A23) 


_L 


_1_ 


oo  R2  (Q1)  *  R(P>Wp(r)  "  awft(r)  * 


CA24) 


where 


wp(r)  * 


[l  +  (1  -  €2)p|2] 


3/2 


Cp' 


[l  +  (1  -  cV2] 
w  (r)  «  ^ 


3/2 


[p"-C2(p"-p'/r)] 

Under  the  conditions  of  practical  interest  mentioned  earlier  Eq.  (A24)  reduces  to 


-L  . 


i_a  _i_  1 


«  R2<Q’>  r(p)  a 


(A25) 


since  Eq.  (A22)  which  now  reads  p(r)  *  r  -  a  *  h  gives  rise  to 


p'(r)  «  1  and  p"(r)  *  0 


(A2K) 


so  that  we  have  wD  *  1  and  w  «  1,  In  Eq.  (A25)  we  see  again  the  conservation  of 
the  relative  curvature  between  the  actual  and  the  RS  geometries.  We  might  men¬ 
tion  that  this  relation  has  sometimes  been  introduced  into  the  literature  without 
due  statement  about  its  condition  of  validity.  We  believe  it  worth  noting  that,  like 
Eq.  (17)  and  the  relation 

R2(Q')  *  -R^P') 

that  follows  from  comparing  Eqs.  (A.  17)  and  (A25),  it  is  of  only  limited  validity. 
The  relation 

RAP')  *  -RJQ')  * - - - 

1  2  l/R(P)  -  1/a 


may  be  re-expressed  in  various  familiar  forms.  Thus,  by  defining  the  lapse  rate 
of  the  refractive  index  byo^  -  -n'  and  noting  that  n  *  1  and  cos  0  8  1  for  all 
practical  purposes  we  have 


m 


where  use  has  been  made  of  Eq.  (A  19).  For  the  special  case  of  n'  ■  const,  the 
above  expressions  for  Rj(P')  reduce  to  those  given  in  Ref.  2.  When  Eq.  (A13) 
is  used  we  have 


R2(Q,)  *  a  •  i  +  an'lcos07n  5  a  *  k<r)  *  ae<r)  ‘ 

In  the  case  of  n'  =  const,  the  approximately  constant  values  of  ag(r)  and  k(r), 
written  as  ae  and  k  and  called  the  equivalent  earth  radius  and  the  earth  radius 
factor,  respectively,  are  of  considerable  utility  in  the  discussion  of  wave  propa¬ 
gation  problems  according  to  the  RSM.  In  particular,  when 

l/4a  *  40  N  unit/km  which  is  customarily  considered  to  correspond  to  the 
so-called  standard  atmosphere,  we  obtain  the  familiar  value  of  ag  *  (4/3)a.  Before 
discussing  the  RS  coordinate  system  employed  in  Ref.  3  we  should  point  out  the 
relation  between  the  elevation  angle  0(r)  in  the  actual  geometry  (see  Figure  A  1(a)) 
and  the  corresponding  angle  called  £(p),  measured  positive  in  the  direction  from 
the  earth's  surface  toward  the  ray  trajectory,  in  the  RS  geometry  (see  Figure 
Al(c)).  We  have  then 


With  i'( r)  given  by  Eq.  (A20)  and  p'(r)  *  1  we  obtain 

tan  6(p)  *  — •  *  sin  6(r) 

^  rn 

For  small  angles,  this  is  equivalent  to  0(p)  *  6(r).  That  is,  under  the  conditions  of 
our  interest  we  may  identify  the  angle  0(p)  with  the  experimentally  measurable 
angle  9(r).  Recalling  the  similar  relation  9(h)  *  9(r)  for  the  EF  system  we  see 
that  we  may  use  the  measured  take-off  or  arrival  angle  as  the  corresponding 
angle  in  either  the  EF  or  RS  system.  In  the  remainder  of  this  subsection  we 
will  summarize  the  RS  transformation  employed  in  Ref.  3  by  slightly  rephrasing 
it  to  make  it  easier  for  us  to  understand.  Consider  a  transformation  which  maps 
a  point  Pwith  polar  coordinates  (r,^)  in  the  real  geometry  [Figure  A  1(a)]  into 
another  point  Pg  with  new  polar  coordinates  (r  ,  ?e>.  Under  general  circumstances 
such  a  transformation  from  (r,  4)  to  (re>  4£)  would  be  represented  by  the  functional 
relations  rg  =  F(r,  <t>)  and  =  G(r,  $).  where  for  convenience  we  introduced 
4  =  »/2  -  ?  and  <t>e  =  t/2  -  Let  us  specialize  the  transformation  to  the  case 
where 

re  =  F(r) 

*e  =  Ffe)  *  * 

Here  F(* )  is  a  function  the  form  of  which  is  yet  to  be  determined.  The  first  equa¬ 
tion  tells  us  that  a  circle  in  the  original  geometry  is  transformed  into  another 
circle  in  the  new  coordinate  system.  In  the  second  equation,  a  stands  for  a  con¬ 
veniently  defined  earth  radiuB,  not  necessarily  equal  to  the  customarily  adopted 
value  of  about  6370  km,  such  as  the  radial  distance  up  to  the  height  of  the  layer 
interface.  What  the  second  equation  means  physically  is  simply  that  the  arc 
length  s  *  a$  and  se  *  F(a)$e  in  the  original  and  transformed  systems,  respectively, 
arc  equal.  If  we  note  the  geometrical  relations 


(A28) 


tan  9  (r  )  = 
e  e 


r  d$  (r  ) 
e  e  e 


and  demand  that  the  angle  $e(re)  =  5e(F(r))  =  9(r)  be  identical  to  0(r),  then 
Eqs.  (A27)  and  (A28)  lead  to  a  differential  equation  for  the  unknown  function  F 
given  by 


rF'(r)F(a)  _ 
aF(r) 


(A29) 


This  has  to  be  solved  under  an  appropriate  boundary  condition.  One  way  to  do  this 
is  to  regard  Eq.  (A29)  as  a  differential  equation  for  a  function  defined  by 
f(r)  *  rF(r).  Its  solution  is  of  the  form 


f(r)  «  (cr)2->  .  (A 30) 

where  c  is  the  integration  constant  and  y,  given  by 


2 


is  another  constant  which  characterizes  the  atmosphere,  as  will  be  seen  below. 
The  constant  c  can  be  given  in  terms  of  y  by  using  the  expression  for  f(a)  obtain¬ 
able  from  Eq.  (A30).  In  this  way  we  can  determine  F  as 

■  nr.  ■  t*7  (*}’  ■  t4t  (*r  • 

For  r  *  a  the  above  equation  gives 

a„  *  F(a)  *  -r-2—  .  (A 32) 

e  1  *  y 

So  far  there  is  nothing  in  the  above  transformation  to  tndicate  that  the  ray  trajec¬ 
tory  in  the  new  system  will  be  a  straight  line.  To  see  that  it  can  be  under  appro¬ 
priate  conditions  let  us  recall  the  ray  equation  for  the  original  geometry: 
rn(r)  cos  ®(r)  *  const.  First,  note  from  Eq.  (A31)  that  we  have 


If  the  refractive  index  is  restricted  to  a  class  of  functions  of  the  form 


n(r)  =  n(a) 


<A33> 


then  the  ray  equation  in  the  original  system  gives  rise  to  the  relation 

r  cos  0  =  const 
e  e 

in  the  transformed  system,  where  use  has  been  made  of  the  relation  0(r)  =  6  (r  ). 
Note  that  the  relation  just  derived  represents  the  equation  for  a  straight  line. 

Thus,  we  see  that  the  transformation  from  (r,d)  to  (r  ,< t  )  that  has  been  employed 
will  describe  an  RS  transformation  for  a  class  of  refractive  index  of  the  form  of 
Eq.  (A33).  A  question  remains:  What  constant  value  of  the  refractive  index  will 
cause  propagation  along  a  straight  line?  To  answer  this  question  we  recognize 
that  a  refractive  index  of  the  form  of  Eq.  (A33)  can  make  practical  sense  only  when 
r  *  a.  In  that  case  we  have  (a/r)"^  *  1  regardless  of  the  value  of  7  and  hence 
n(r)  «  n(a).  That  is,  choosing  almost  any  value  of  n(r)  for  which  r  *  a  will  make 
no  difference.  In  Ref.  3  the  authors  use  the  value  at  the  height  of  the  layer  inter¬ 
face  to  calculate  the  time  delay,  which  will  be  perfectly  acceptable. 

The  significance  of  the  constant  7  that  characterizes  the  atmosphere  may  be 
found  in  the  following  way.  For  an  atmosphere  with  constant  gradient  of  the 
refractive  index  we  have  for  r  -  a  *  h  «  a. 


a.C 

n(r)  *  n(a)  +  n'(h)  «  n(a)(l  -«^h)  *  n(a)  (-p-j 


Comparing  this  with  Eq.  (A33)  we  fi.it  that 


Having  examined  the  transformation  employed  in  Ref.  3  we  may  observe  that 
it  represents  a  subset  of  the  more  general  RS  transformation  discussed  in 
Section  A2.  3.  We  say  this  because  the  transformation  is  valid  only  for  the  case 
where  the  transformed  surface  of  the  earth  has  a  constant  radius  of  curvature, 
whereas  the  radius  RgfQ')  in  Section  A2.  3  need  not  be  constant.  Only  when  the 
ray  trajectory  remains  almost  horizontal  and  close  enough  to  the  earth's  surface, 
and  when  the  refractive  index  has  a  constant  gradient,  does  our  general  RS  system 
reduce  to  the  one  used  in  Ref.  3.  Since  we  will  be  concerned  exclusively  with  the 
circumstances  in  which  the  conditions  just  mentioned  prevail,  the  RS  transforma¬ 
tion  in  Ref.  3  will  be  seen  to  be  an  appropriate  one  to  be  employed  for  the  analysis 


Appendix  B 

Comparison  of  the  EF  and  RS  Descriptions 

1*1.  THE  EFM 

In  this  method  the  earth's  surface  is  flat  and  the  ray  trajectory  is  curved, 
for  an  atmospheric  layer  with  a  constant  refractive-index  gradient,  n'  =  -  = 

const,  the  ray  trajectory  is  given,  under  the  assumption  that  r  =  a  and  cos  0*1, 
by  a  circular  arc  whose  radius  is 

R  a  - - - •  (Bl) 

oC  -  1/a 

The  elevation  angle  .0(h)  at  an  arbitrary  point  on  the  ray  trajectory  in  this 
system  is  related  to  the  elevation  angle  9(r)  in  the  actual  system  (see  A2.2)  by 

tan  0(h)  =  —  tan  0(r)  I 

a  . . 

|r=a+h 


Therefore,  under  the  conditions  of  our  interest  we  may  set 


B2.  THE  RSM 


Here  the  ray  trajectory  is  represented  by  a  collection  of  straight  line  segments 
that  may  be  rectified  into  a  single  straight  line,  while  the  earth's  surface,  layer 
interface,  etc, ,  and  other  lines  are  transformed  into  curved  lines.  According  to 
the  particular  version  of  the  RS  system  employed  by  the  Signatron  workers  the 
earth's  surface  (radius  a)  and  the  layer  interface  (radius  r^  =  a  +  h^)  in  the  plane 
of  the  ray  trajectory  are  transformed,  respectively,  (see  Eq.  (A3 1  > J  into  circles 
of  radii 


a  =  F(a)  =  — —  and 

e  1  -  -y 


reD 


F(rD)  - 


(B3) 


The  elevation  angle  0  (r  ),  at  an  arbitrary  point  on  the  ray  trajectory  in  the 
RS  system,  is  exactly  equal  to  the  elevation  angle  0(r)  at  the  corresponding  point 
in  the  actual  system  (see  below  Eq.  (A28)(  : 

S  (r  )  =  0(r>  .  (B4> 

e  e 


B3.  COMPARISON  OK  RESll.TS 

Vte  can  now  proceed  to  examine  an  event  occurring  in  the  actual  system 
(Figure  Bl(a))  as  seen  in  the  EF  and  Kb  systems  (Figures  Bl(b)  and  Bl(c)). 

For  the  EF  system  we  have  from  Figure  Bl(b)  and  Eq.  (B2) 

Lj  *  Retain  -  sin  0j) 

V  *  Hf%  *  V  •  (B5) 


On  the  other  hand,  from  Figure  R 1  (c  >  we  have 


I  '  *  F(r  1(0 
1  ,rl)  '  1) 


-  V 


X  i  =  F(rjj)  sin  0^  -  F(r.p)  sin  6-j. 

where  r  j.  *  a  +  h^,.  Bv  using  the  relation  F(r>  =  1  "  fl-p  we  ma.v  writ« 

X  *  Ffr^Hsin  -  sin  11  j.)  *  d-j,  sin  *  «  here  the  first  term  on  the  right  hand  suit 

equals  Xj'.  For  h[:)/a  «  1  we  have  K(r(y  -  Ka)(a  r^)"''1  *  F(a>  -  a{  so  that  w. 
may  write 


9. 


» 


9, 

S-: 


i.-4 


‘i 


■  i 


a. 
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Rome  Air  Development  Center 

*AK  pla ns  and  extent*  research,  development,  test  and 
selected  acquisition  programs  in  s upport  of  Command,  Co, 'trial 
Communications  and  Intelligence  (C*l)  activities.  Technical 
4appo«  within  areas  oi  technical  competence 
as  provided  to  ESV  Program  Offices  (POs)  and  other  BSD 
elements.  The  principal  technical  mission  areas  axe 
communication*,  electromagnetic  guidance  and  c ontrol,  6ur- 

wetra!^e  °*  ***o*pace  object s,  intelligence  data 

collection  and  handling,  information  system  technology, 
ionospheric  propagation,  4 olid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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